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The Mathematics of the Physical Properties of Crystals 

By WALTER L. BOND 

SECTION 1 
Introduction 

THE use of crystals as oscillating elements and as light valves in electric 
circuits has given the mathematics of crystalline media an engineering 
importance. Soon after the first simple quartz oscillators were made it was 
noticed that some ways of cutting the block from the natural crystal gave 
lower temperature coefficients of frequency than did other ways. This led 
to studies of the change of elastic modulii with direction and temperature 
and finally to the discovery that there are directions in quartz for which the 
shear modulus does not change with temperature. 

Such computations are rather involved, and there is, in the English 
language, no general reference book on these new problems. The existing 
works were evidently not written with the idea in mind that anyone would 
ever actually do much numerical work with directional properties of crystals, 
since the methods used are not the best suited to this. The matrix algebra 
has the advantages of a symbolic algebra and is also, through the concept of 
matrix multiplication, a scheme for computing numerical results. 

As the problem of temperature coefficients of frequency involves the 
temperature coefficient of expansion, the temperature coefficient of density 
and the temperature coefficient of elastic modulii, these problems must be 
put into the language of matrix algebra so that they will fit into the general 
structure being built for more difficult problems. For this reason, after an 
introduction to the idea of linear vector functions, through consideration 
of the relation between the electric field and the induction in a crystal, and 
a hasty sketch of symmetry types found in crystals, we proceed to the 
consideration of stress and strain and their relations to each other. 

Following these, we take up piezo electricity and the converse piezo 
electric effects; these are important as they tell us the ways a crystal may be 
driven. We have not seen anywhere a general proof that the modulii of 
the converse effect are the same numbers as the modulii of the direct effect 
— to the first order of small quantities, though Lippman predicted the 
converse effect and demonstrated its magnitude to be about this; he ap- 
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parently didn't consider the general case, of six simultaneous stress com- 
ponents, six strain components, three electric field components and three 
induction components. The fact that the mentioned relation is true only 
to the first order of small quantities seems to have escaped the attention of 
some experimenters who have sought to show non-linearity of the piezo- 
electric effect by demonstrating non-linearity in the converse effect. 

As a basis for light valve problems, we handle the propagation of light 
through crystals, then the electro optic effect and the piezo optic effect. 

SECTION 2 
A Linear Vector Function 

For almost every physical constant of an isotropic medium a crystalline 
medium has several constants. For instance, a piece of glass has a co- 
efficient of thermal expansion but a crystal has many coefficients of thermal 
expansion, the coefficient depending on direction. It might be thought 
that there were no necessary relations between the coefficients in different 
directions but there are necessary relations. 

As an example of the simplifying relations between the values of physical 
constants in different directions let us consider a crystal in an electric field. 
Measurements show that the dielectric constant varies with direction in a 
crystal. If the field is not in the direction of greatest dielectric constant, 
the displacement current might veer over a little, much as a nail tries to 
follow the grain of the wood. We shall assume that for any electric field 
vector E there corresponds an electric induction vector D which may not 
coincide with E. Also we assume that the magnitude of D is proportional 
to E, that is, if E results in D, then nE results in nD. Lastly, we assume 
that if Ei results in D h E% in L\ . . . and E m results in D m , then Ei + E 2 \ . . 
E m results in Z>, + A + . • • D m . If these assumptions hold, then as any 
arbitrary field E can be expressed as the vector sum of its three components 
Ei, E2, E 3 along three arbitrary .unit vectors i, j, k, the induction vector 
resulting from E can be computed from the induction vectors resulting from 
Ei, E2, and E 3 . For, let Ei result in D u i + D*j + D 3 ik, E 2 result in D n i + 
D-aj + Dmk and E 3 result in D u i + Atf + D^k, then E x + £2 + E 3 = E 
results in the induction vector: 

D = (EiDu + Ed> n + EMi 

+ (.E1A1 + EsDa + EM)j (2.1) 

-I- (EiDn + E»Du + -EsAa)* 

It is seen then, that not more than 9 constants are needed to describe 
the dielectric properties of a crystal. The energy required to establish the 



MATHEMATICS OF PHYSICAL PROPERTIES OF CRYSTALS 3 

electric field is half the product of the component of the induction in the 
direction of the field and the electric field. This is, therefore: 

2W = E\D n + E\D n + ^sAk + ^^(As + A2) + ^(Ai + #13) + 

-Ei-E 2 (A 2 + Ai) 

Considering then a condenser made from a unit cube of crystal, the charge 
is D and the energy content is IV. If there is no leakage loss, the charge 

that can be drawn from the condenser is D = — - . Whence A = — f = 

ah dE\ 

A1-E1 + KA2 + Ai) E 2 + \ (As + Ai) E 3 . If, therefore, the induction 

is derivable from a potential, A2 = \ (A2 + Ai) or D a = Ai- Similarly 

A3 = Ai and A3 — A2. By a proper choice of axes the remaining six 

D's can be reduced to three. In the case of isotropic dielectrics Ai = A2 = 

A3 and 4ttAi corresponds to k, the dielectric constant. 

SECTION 3 

The Symmetry of Crystals 

If a crystal has certain sorts of symmetry the number of constants re- 
quired to describe each property is materially reduced. For this reason 
we now turn our attention to a study of symmetry. 

In general, plotting a vector property of the medium for a crystal gives a 
complicated surface which we shall call a property surface. Each property 
surface of a homogeneous isotropic medium is a sphere. 

Because of the orderly arrangement of matter in a crystal, the property 
surfaces of crystalline media are commonly symmetrical. If a casting of a 
property surface were made it might fit into its mold in several positions. 
A property surface for quartz for example, if lifted from its mold and rotated 
through a third of a turn about the proper axis, would fit back into the mold. 
That is, quartz has a three fold axis. The natural requirement that mole- 
cules be laid down in a way economical of space limits the kinds of symmetry 
possible for crystals to axes of two fold (binary) symmetry, of three fold 
(trigonal), of four fold and of six fold symmetry, planes of reflection symme- 
try and combinations of axis-reflection symmetry, besides a simple sym- 
metry through a center. From these elements it is possible to divide all 
possible property surfaces into 32 classes. No other classes built from these 
elements could be self-consistent. 

A diagram study will prove this point. On a sphere let us mark axes of 
two fold symmetry by means of a solid boat shaped figure, three fold with a 
solid triangle, four fold with a square, six fold with a hexagon, planes of 
symmetry with a solid line (great circle) and combination axis reflection, 
by means of similar hollow figures. Finally, we shall project the sphere 
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FIG. I 
ASYMMETRIC 



FIG. 2 
PINACOIDAL 



MONOCLINIC SYSTEM 



FIG. 3 
SPHENOIDAL 





ORTHORHOMBIC SYSTEM 



H-+-T* 



FIG. 6 
BISPHENOIDAL 




FIG. 7 

PYRAMIDAL 



FIG. 8 
BIPYRAMIOAL 



TETRAGONAL SYSTEM 




FIG. 9 
BISPHENOIDAL 



FIG.IO 
PYRAMIDAL 



FIG. II 
SCALENOHEDRAL 



FIG. 12 
TRAPEZOHEDRAL 




FIG. 13 
Bl PYRAMIDAL 




FIG. 14 

DITETRAGONAL 

PYRAMIDAL 




FIG. 15 
DITETRAGONAL 
BIPYRAMIDAL 



Figs. 1-15 — Crystal classes. 
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TRIGONAL SYSTEM 




FIG. 20 
DITRIGONAL 
PYRAMIDAL 



FIG. 21 

DITRIGONAL 

SCALENOHEDRAL 



FIG. 22 
DITRIGONAL 
BIPYRAMIDAL 



FIG. 23 
PYRAMIDAL 




HEXAGONAL SYSTEM 
\ o/ ! \x> 



FIG. 24 
TRAPEZOHEDRAL 





FIG. 25 
BIPYRAMIDAL 



FIG. 26 
DIHEXAGONAL PYRAMIDAL 



FIG. 27 
DIHEXAGONAL BIPYRAMIDAL 



/x 4 "jO^O \ 

f — 8-M — • 



FIG. 28 

TETRAHEDRAL- 

PENTAGONAL- 

DODECAHEDRAL 




CUBIC SYSTEM 

FIG. 29 

PENTAGONAL 

COSITETRAHEDRAL 





FIG. 30 

DYAKIS 

DODECAHEDRAL 



FIG. 31 
HEXAKIS-TETRAHEDRAL 



FIG. 32 
HEXAKIS- OCTAHEDRAL 



Figs. 16-32 — Crystal classes. 
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and markings onto a plane through the center. Figures 1 to 32 is a set of 
such diagrams. Fig. 23 for instance shows a six fold axis. Fig. 1 represents 
a medium with no symmetry whatever. The cross represents a typical 
vector property, the vector piercing the sphere above the projection sheet. 
If the vector pierced below the sheet it would be marked with a circle. The 
dashed circle of Fig. 23 indicates the boundary of the sphere without im- 
plying it to be a plane of symmetry. The presence of six fold symmetry 
requires the typical vector to be shown in six places. If an axis of two fold 
symmetry is added at right angles to the six fold axis, it must appear six 
times and the typical vector must now appear twelve times, six times above 
and six times below the projection sheet. Continuing in this way we shall 
find the self-consistent classes of symmetry to be the 32 shown in the dia- 
grams. Often the symmetry of a crystal class is expressed by means of a 
formula. A center of symmetry is symbolized by the letter C, a binary 
axis by A 2 , a trigonal axis by A 3 , a ternary axis by A it a six fold axis by A 6 , 
a plane of reflection by P, and a combination rotation reflection by the 
combination symbols A> 4 or AV In this way the symmetry formula of 
quartz for example, is 3.4 2 - .4 3 . 

SECTION 4 

Matrix Algebra 

In the solution of problems of crystal physics we are involved in the 

handling of many sets of linear simultaneous equations. As the matrix 

algebra lessens the work involved in handling sets of linear simultaneous 

equations we turn now to a study of matrix algebra. 

Several independent variables Xt, x 2 . . . x n are linearly related to several 
other independent variables yi, y 2 . . . y m as 

y\ — flllS-'l 4" 012*2 + • • • dlnXn 
y 2 = 021*1 + • • • 



= dmlXl + • • • dmnXn 



or briefly 



yi 



= £<%** 1-1,2 ••• m (4.1) 



j-i 



In most all such equations as (4.1) the variable to be summed over appears 
twice in the subscripts of one side. As a convention we agree to omit the 
summation sign and sum wherever subscripts are repeated. 
Thus: yi = aax,- is to be summed over/ 
again, if x, = bj k z k the 2's being a third set of variables we have: 

y t = aijbjkZk to be summed over j and k. 
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We can think of this as a special multiplication of hyperquantities a, b 
and z. If we define 

Cik = dijbjk (4.2) 

we may go from the y's to the z's directly thru y,- = e«8*. We can now 
consider the "table" 

'flnOi2 • • • a n ^ 
(hi • • • 
\a m i . . . c mn 
as being the quantity a, and the table 

/buba • • • • 



\b* 



as the quantity b. 

These "tables" are called matrices. 

Going to eq. (4.2) we see that the quantity c is to be a "table," the typical 
element c,-,- of which is to be gotten by multiplying the ith row of a by the 
yth column of b, term by term thus: 

Cij = anbij + Oahj + • ■ • 

After a little practice it becomes almost automatic to form the yth term 
of the product of two matrices by letting the index finger of the left hand 
follow across the ith row of the left matrix while the right index finger 
follows down the ,/th column of the right matrix. The fingers step along in 
synchronism and at each pause the quantities under the two fingers are 
multiplied and the product added algebraically to the accumulated sum. 

The algebra of these special multiplications is not commutable, i.e. 
ab 7± ba. 

Eq. (4.1) can be considered as a special case of eq. (4.2), in which the 
matrices x and y have one column only. In this manner a vector with 

components x\ x% .r 3 can be considered as the matrix 

If eq. (4.1) has the same number of .v's as y's we may solve (by means of 
determinants) for the x's in terms of the y's. We would then get a new set 
of equations 

*< = (hi y,- 



8 BELL SYSTEM TECHNICAL JOURNAL 

The significance of the a" 1 is that the matrix product of a and a~ is a matrix 
with ones on the major diagonal and all other terms zero. Whenever the 
product of two square matrices gives such a matrix (known as the idem- 
factor, J) they are said to be reciprocal. Only square matrices have 
reciprocals. Multiplying any matrix by the idemfactor leaves the matrix 
unchanged. We might consider, as part of our mathematical short hand, 
that eq. (4.1) was solved for x by multiplying through by a~ , as 

a~ l y = aT l ax = Ix = x. 

We must remember that the order must not be disturbed as the quantities 
are not commutable, and that only square matrices have reciprocals. 

The major diagonal of a square matrix is the set of terms running diagon- 
ally from the upper left to the lower right. 

A symmetrical matrix has any term M,-,- = M ',■; 

An anti-symmetric or skew symmetric matrix has any term M ,-,• = — Ma 
for i 9^ j. 

Rotation Theory 

The matrix algebra can be used to express a vector as a function of another 
vector, that is to handle such relations as exist between E and P of section 2. 

There is another important aspect of matrix multiplication, that of trans- 
forming a function from one set of axes to another. Let us assume that the 
new set of unit axes, x[ x' 2 and x' 3 are merely the old ones rotated through 
angle <t> about some axis A which is a unit vector passing through the origin. 
From Fig. 33 we see that in the expression: 

Xi = 0-11*1 + 021*2 + 031*3 

the aa's are the cosines of the angles between xi and the three quantities 
x'j. Conversely they are the cosines of the angles between the ay's and %i. 
Consequently, if the primed unit vectors are given in terms of the unprimed 
ones by the three equations 



x t = aijXj 



then the unprimed x*s are given in terms of the primed ones by the ex- 
pression: 

Xj "~~ ai\X% 

This reversible relationship is well depicted by the table: 





-1-1 


Xi 


*3 


Xl 


flu 


«12 


013 


X2 


021 


022 


023 


1 


«31 


032 


fl33 
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In this direction cosine table we can "look up" the components of any unit 
vector in terms of the other system. 

The matrix a a is merely the matrix a,-,- with rows and columns inter- 
changed, an is called the conjugate of a,,-. We shall denote the conjugate 
of any matrix M by M c . 

Obviously V is the vector sum of the 3 components (on the new system) 
of each of its 3 components on the old system. 

ZaTiVi + aTiV* + al}V»\ 
(V) aew - [a* V + 




Fig. 33 — The direction cosines of Xi on X'i X' 2 XV 

If the expression giving the components of V on the new system is de- 
noted by V we may write 

V = aV 

conversely V = a -1 V 

Since Xi is of unit length, the sum of the squares of its three components 
(on the primed system) is unity. 

That is On + o 12 + a i3 = 1 

(4.3) 



similarly 
and 



2 2 

#21 + ^22 + #23 = 1 

^31 + 032 + 033 = 1 



Now a c can be considered as a rotation similar to the rotation a. Con- 
sequently their product aa c is a similar rotation. Let us consider this 
product. 



~ ! (4.4) 
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The squares of its terms must sum to zero, row by row as in (4.3) 

(flllfll2fll3\ /flllfl2l03l\ /l • A 

021022^23 J I 012022032 1 = I • 1 • I 
031032033/ \oi3O23033/ \. . 1/ 

Because of the relations a n + o? 2 + flis = 1, etc., we see that the terms of 
the third matrix are zero for all terms not on the major diagonal. There- 
fore, aa c is an idemfactor and the reciprocal matrix of o is the same as its 
conjugate matrix. 

o c = o 

Also x[ is of unit length, and the sum of the squares of its components on 
the unprimed system is unity. Thus we find: 

2.2,2 1 

On + a 2 i + a 3 i = 1 

012 + 022 + O32 = 1 (4.5) 

2 I 2 1 2 _ a 

013 T 023 "T O33 — 1 

We now introduce from vector analysis the concept of the scalar product. 
The scalar product of two vectors u and v is u c v. It is the product of the 
lengths of the two vectors and the cosine of the angle between them. 

If we take the scalar product of Xi and x 2 as expressed in the primed 
system we have, since they are mutually perpendicular: 

/fl2l\ 
(flll, O12, Ois) I O22 I = Oll021 + 01202? + Oi3fl23 — 

Similarly multiplying x 2 and x 3 scalarly, and x 3 and .Ti we find: 

Allfl21 + 012022 + Ol3fl23 = 

02lfl31 + 0220-32 + 023033 = (4.5) 

fl3lOll + 032012 + 033013 = 

If we multiply #i and x 2 etc. as expressed on the unprimed system we get 
the relations: 

OllOl2 + 021022 + 03lfl32 = 

fll2fll3 + 022fl 2 3 + 032033 = (4.7) 

013011 + 023021 + 033031 = 

The vector product of two vectors u and v requires the defining of a 
special matrix, the cross matrix. 

(0 — «3 «2\ 
u» -Mil (4.8) 
— «2 Wl 0/ 
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We note that this is formed by writing zeros on the major diagonal, then 
going back, from the lower right corner writing u h «2 and u 3 around the edges. 
We then make the lower left term negative, then operate on the opposite 
side of the major diagonal so as to make the matrix skew symmetric. 

The reciprocal of any matrix m is 

^=«' + ' (4.9) 

where M a is the ji minor of | m \ . 

The cross matrix has no reciprocal as for it (4.9) becomes indeterminate. 

Since the vector product of two vectors u and v is another vector per- 
pendicular to both u and v and of a length uv sin (uv) we may write, on the 
primed system 

*i X ff2 = xz in the form 

/ — an 0iA AfeA / — 013022 + an(hz\ 

I 013 —an I I 022 1 = I #13021 — flll023 I = 

\ — fll2 flu / \fl23/ \ — Al2fl21 + 011022/ 

Matrices including vectors are equal only when their corresponding terms 
are equal. Hence, we get the relations 

031 = Al2023 — Al3fl22 

flB2 = Al3fl21 — Allfl23 • (4.10) 

fl33 = Allfl22 — Al2021 

Similarly we get the relations: 

All = A22fl33 — A23fl32 

012 — A23fl31 — A21033 

013 — A2lfl32 — fl22fl31 f. **\ 

(4.11) 

021 = A32013 — 012033 
fl22 = fl33flll — A3lfll3 
fl23 = fl3lOl2 — fl32flll 

The 21 relations between the fl,/s allow us to complete the matrix given 
four terms. 

Several Useful Matrix Relations 

/B/d xA 

The del operator is the pseudo vector V = I d/d .r 2 J (4.12) 

\d/d xj 
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It transforms upon a rotation of axes as does an ordinary vector: 

V' = aV (4.13) 

grad u = V u c , a matrix (4.14) 

div u = V c u, a scalar (4.15) 

curl u = V X «, a matrix (4.16) 

grad radius vector = Vp = I, the idemfactor (4.17) 

(abc...)- l = ...c-W (4.18) 

(o&c . . .) c c = ... c c b c a e (4.19) 

'a n ..A" 1 /1/fln . 

022 ...\ = 1/022 ... , , 

33 1/033...' 



,« , . • \-i (matrix) /,•««% 

(Scalar tunes matrix) = r-^— (4.21) 

x scalar 

SECTION 5 

The Geometry of Rotations 

As a first application of the matrix algebra let us compute the a matrix 
for a few general rotations. Although we can consider a general rotation 
as one of angle <f> about the unit vector* A, it is easier to consider a general 
rotation as three successive rotations about coordinate axes. 

A study of Fig. 34 shows that for a counterclockwise rotation 4> about %\, 
the new components of a vector V are: 

V[= V x 

y'i = Vi cos </> + Vz sin 
V 3 = —V 2 sin 4- V 3 cos <f> 

whence V = aV where 

/I \ 

c = cos sin 1 (5.1) 

\0 — sin cos 0/ 

* A general rotation of amount about the unit axis A is given by 
a = AA e + (I - AA C ) Cos <f> + Sin A 
See Vector Analysis (Gibbs Wilson, Yale Press) pp. 338. 
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Similarly, for a counterclockwise rotation 4> about %% we have 

(cos — sin <f>\ 
1 1 (5.2) 

sin (f> cos 0/ 

and for a counterclockwise rotation about X3: 

(cos </> sin 4> 0\ 

— sin cos 4> I (5.3) 

1/ 
x 3 




Fig. 34 — The relationship between the components of a vector on one coordinate 
system and on another. 

In the appendix we give the special transformations corresponding to the 
symmetry operations of the 32 crystal classes. If we have three successive 
rotations: 

x' = ax 



the resultant rotation is 

or 
where 



x = ax 
x'" = a"x" 



%'" = a" a' ax 
x'" = Rx . . 
R = a" a' a. . 



(5.4) 
(5.5) 
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The I.R.E. Orientation Angles and the I.R.E. Matrix 

The Institute of Radio Engineers has proposed that, for quartz crystals, 
all orientations be given in terms of three rotations <f>, 6, ty about x 3 , x% 
and x 3 respectively, starting with the plate length along X\ width along 
X2 and thickness along x 3 . (Here x 3 is the z or optic axis, X\ is the electric 
axis.) 
Whence, here: 

cos <f> sin ()\ 
— sin <f> cos <f> J 
.0 1/ 

and carrying out the two matrix multiplications: 




Xi X2 x s 

cos <f> cos cos if/ sin <£ cos cos ^ — sin cos \p 

— sin </> sin ^ +cos<£sin^ 

— cos <f> cos sin ^ — sin <f> cos sin \p 

— sin <f> cos \j/ +cos <f> cos ^ sin d sin \p 



cos ((> sin 6 



sin (j> sin 



cos 



Xi 



Xi 



in 

Xz 



...(5.6) 



If we denote the unit vectors along the length, width and thickness as P\Pi 
and P 3 respectively we have as a matrix defining the plate: 

P = Rx (5.7) 

The I.R.E. orientation system is useful to the designer of crystal plates 
because his problem is to choose such values of </>, 0, ^ as to give the plate 
certain physical properties along its length, width and thickness. The man 
who cuts the plate has a different problem, that of moving the crystal (and 
hence the xi x^ x 3 axes) about a fixed saw so that the plate cut parallel to the 
saw blade is what the designer ordered. 

Let us consider such a system as shown in Figs. 37, 38 and 39. In Fig. 
38 the crystal stands with its optic axis along P 3 , its + electric axis (for 
right hand quartz) along P. Since the shop man considers clockwise rota- 
tion as positive we now rotate the crystal through angle U 3 about P 3 clock- 
wise, we then turn the crystal through angle U 2 clockwise about Pi, and 
finally, after cutting out a slab of required thickness, we turn it clockwise 
through angle L\ about P 3 to cut its length and width. 

On the plate axes Pi the crystallographic axes X\ x% x 3 are now given by 

x = rP (5.8) 
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*3 




Fig. 35 — The initial position (0, 0, 0) for the I.R.E. direction angles. 




Fig. 36 — The final position (<!>, 6, *) for the I.R.E. direction angles. 



. 



16 

where 
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I cos U\ sin f/i 0\ /l 6 6 \ / cos U» sin U 3 0> 

r = ( —sin U\ cos £/ x ) 10 cos £/ 2 sin C/ 2 I I -sin U 3 cos U s 

\ 1/ \0-sin U 2 cosUj \ 1/ 



or 



cos Z7i cos £/ 3 . cos Ui sin t/ 3 sin U\ sin £/ 2 

— sin U\ cos £/ 2 sin Z7 3 +sin U\ cos C/ 2 cos U 3 

— sin t/i cos U 3 — sin £/i sin U 3 cos £/i sin C/ 2 

— cos Ui cos C/ 2 sin U 3 +cos 27i cos Z7 2 cos U 3 



sin t/ 2 sin U 3 



— sin Z7 2 cos I/3 



cos U 2 



(5.9) 



■ 



P2-» 





P,-* 



Fig. 37 



Fig. 38 



Fig. 37— The (0. 0, 0) position of a shop system of direction angles. 
Fig. 38— The second position of a shop system. 

From (5.8) we see that P = rx and hence, if this is to be the same plate 
the designer specified by P = Rx we must have R = r whence we may 
equate the terms of (5.6) and (5.9) to get the relations 
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cos Z7 2 = cos or' #2 = ± 

tan Ui = cot yp 

tan Us = cot <f> 



(5.10) 



or 



U t = $ -90 + mr 

U 2 = d 

Ut = <f> -90 + mr 



P 9 < 



PIEZOID 




Fig. 39 — Cutting the slab and trimming it to the piezoid boundaries. 



18 ■■•.' BELL SYSTEM TECHNICAL JOURNAL 

Where n is any integer positive or negative, including zero. If we take 

V t - ^ - 90 

U 2 =d (5.11) 

U 3 = + 90 

The matrices are consistent term by term. 

SECTION 6 
Crystalline Dielectrics 

As a first application of the matrix algebra considered as a linear vector 
function let us reconsider the problem of the crystal in an electric field. 

The relations of chapter II, equation (1) can be written in the abbreviated 
form: 

D = DijE where D Tt - D tr 

in accordance with the system of abbreviations adopted in the appendix. 
If we put 

4irD Ta = k Ta 



equation (1) can be written 



D = j-K.E (6.1) 

47T 



In order to investigate the effects of crystal symmetry in determining 
the least number of dielectric constants that are required for a given class of 
symmetry it is desirable to find the electric induction D for any system of 
axes. Suppose that we choose a system for axes *i , *2 , * 3 related to *i , 
X2 , %z through the relations: 

Xi = 0n*l + 012*2 + 013*3 

X 2 = 021*1 + <*22*2 + 023*3 (6-2) 

*3 = 031*1 + 032*2 + 033*3 

where a n is the cosine of the angle between X\ and *i , 012 is the cosine of the 
angle between *i and x 2 etc. 

Equation (6.2) can be abbreviated to 

*' = ax 

where is the matrix 
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(flu Ai2 flw\ 
On 022 A23 1 
fl31 A32 033/ 



It is shown in the sec. 4 that any vector V = I V% J can be written on the 




new system of axes as V where V = aV, conversely V = a~ V; a~ is the 
matrix reciprocal to a. Since the induction D and the electric field E are 
simple vector functions they transform as the vector V, that is: 

D' = aD (6.3) 



But by (6.1) 



whence: 



or 



if 



E' = aE (6.4) 



D = 4-kE 

47T 



aD = — aka c aE 
Air 



V =-r- k ' E ' (65) 



k' = aka c , (6.6) 



We see that the form of (6.5) is the same as that of (6.1) for any set of 
axes if (6.6) is used to define the new dielectric matrix k. 

To apply this relation (6.6) to a particular crystal let us consider a tetrag- 
onal crystal (which has its properties unchanged by a rotation of 90° 
about a four fold axis) . Let us choose the four fold axis as Xa and then rotate 
the axis 90° about x 3 . In this case 




I and the reciprocal matrix a J = [ — 1 



10^ 
•100 
1, 



whence equation (6.6) becomes: 



/ 1 0\ (kn 



k a k n \ /0-10\ / fa-kn V 

„*«&») (l 00) = (-£i2 ku-k n 

1/ \kn kn kj \0 1/ \ kn -k n k&, 
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But because of the symmetry k is unchanged by this transformation, 
hence 

k = k' 

Two matrices can be equal only if corresponding terms are equal, hence 

£ll = h 2 , k\2 = — hi = 0, #23 = — &31 and £31 = ^28 

whence 

£23 = £31 = 0. 

We are left then, with the dielectric constant matrix for the tetragonal 
bisphenoidal class: 

/An \ 
k = ( ku.0 ) 

\o kj 

Applying other transformations possible for tetragonal crystals gives no 
further simplification. 

If we go through all the symmetry transformations possible for the 32 
classes we find that cubic crystals require but one dielectric constant, hex- 
agonal, trigonal and tetragonal crystals require two constants, orthorhombic 
monoclinic and triclinic crystals require 3. 

As the triclinic class has no fixed axes or planes of symmetry the reduction 
of its 6 constants to 3 is not so obvious. It may be seen by expanding into 
ordinary xyz coordinates, that p c kp = 1 is the equation of an ellipsoid, (p 
is the radius vector) where the six k's are the coefficients of x , y , z , yz, 
zx and xy respectively. If we choose the coordinate axes along the axes of 
the ellipsoid the yz, zx and xy terms drop out and only three k's are needed. 
With triclinic crystals then, if we determine the axes of the ellipsoid, then 
choose the coordinate axes along them, only three dielectric constants are 
needed to completely specify the polarization in terms of the electric field. 
The determination of the ellipsoid axes must be made experimentally as 
there are no symmetry elements to guide us. It is possible to compute the 
positions of the axes from the 6 k's by solving a cubic equation. 

The values of the k's depend on the frequency of the applied field. In 
crystals of low symmetry the ellipsoid axes for different frequencies do not 
necessarily coincide. 

Another vector quantity of interest is the polarization, P = D — — E. 
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P=^(k-I)E (6.7) 



SECTION 7 

Quadratic Forms 

Often the elements of a matrix are themselves functions of other quanti- 
ties. In order to relate the elements of one matrix with those of another by 
means of a matrix multiplication, we may make a single column matrix of 
each of them. We then wish to know how a transformation of axes changes 
the elements of this single column matrix. Consider a symmetrical matrix 
b that relates two vectors u and v: 

u = bv. 

A transformation of axes, a, changes u and v to u' and v'. Multiplying u = 
bv through by the prefactor a we have 

au = abv. 

We now replace au by its equivalent a~ v' whence: 

u' = aba~W 

so that 

u' = 6V 

if we define b' as 

b' = aba~ 



(7.1) 



To be in accord with common usage we now rearrange b according to the 
arbitrary scheme: 



M 




M 


622 




B 2 


633 


= 


B 3 


^23 






631 






6l2 




kJ 



We wish to know what operation to perform on B to get B' corresponding 
to b'. If we expand b' = aba~ l it is easily seen that b' = aB where 
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2 

flu 

2 

(hi 

2 

o 3 i 



2 

2 
022 

2 

a 32 



2 
Ol3 

2 

(ha 

2 

033 



fl21 031 022 ^32 023 O33 



O31 On 032 O12 O33 C13 



On 021 Ol2 022 O13 023 



2fli2fll3 20i3 0ll 2flnOi2 

2022 023 2023 021 2021022 

2032 033 2033 031 2031032 
022 033 021 O33 O22 O31 

+ 023 032 +O23O31 + 021032 

O12 O33 O13 O31 On O32 

+ 013 032 +O11O33 +O12O31 

012 023 013 021 O11O22 

+ 013 022 + On 023 +O12O2I 



■ -(7.2) 



Because we shall often need to form the a matrix from the a matrix we 
need an easily remembered mechanism for doing so. We notice that those 
are four kinds of terms in the a matrix and that the four kinds can be sepa- 
rated from each other by two center lines, one horizontal, one vertical. 
This gives us four squares of nine terms each and we can correlate each term 
of any square to a term of the a matrix by means of its position in the square. 
The terms of the upper left square are the squares of the corresponding terms 
of the a matrix. To form any term of the lower left square we cover the 
corresponding term of the a matrix with our finger and multiply the visible 
terms of that column. To form any term of the upper right square we cover 
the corresponding a term and write down double the product of the visible 
terms of that row. To form any term of the lower right square we find the 
corresponding o term, strike out that row and column and write down the 
sum of the remaining cross products. A study of the following diagram 
will help to remember these rules. 



Terms are squares 

of corresponding 

a terms 


Omission products 
doubled 


Omission products 


Sum of omission 
cross products 



Fig. 40 



SECTION 8 

Crystal Elasticity 

Stress 

Consider a point P in a medium acted on by forces. If a small area is 
chosen about P the medium on one side of the area exerts a force on the 
medium on the other side. The force will depend on the size of the area and 
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on the direction of its normal n. We shall choose a triangular area ds such 
that an arbitrarily chosen set of mutually perpendicular unit axes Xi , xz , 
Xi pass through the vertices of the triangle. Let us consider the conditions 
of equilibrium of the tetrahedral element of volume so formed. The areas 
normal to Xi , x% , x 3 are dsi , ds 2 , ds 3 , respectively, and the forces per unit 
area acting through these faces are: 





Any body forces (such as gravity) depend on a higher order of smallness 
(that is on the volume rather than on the area) and hence are negligible. 
Whence for equilibrium: 



But 



Fds = F x dsi + F 2 ds 2 + F 3 ds 3 
dsi = ttids, ds 2 = n 2 ds and ds 3 = n 3 ds 



where I w 2 I is the normal to the area ds. Whence we may write: F = fn 

w 

where/ is the matrix 

/fn /12 /i3\ 

I /21 /22 /23 1 

\fnj 32/33/ 

For the body to be in rotational equilibrium the tangential forces must 
balance, hence / M = /21 ,/i 3 = /31 and/ 23 = / 3 2 • 

Transformation of Axes 

A change of axes that transforms vectors through F' = aF changes F = 
fn to aF = afa~ an so that if /' = afaT 1 then F' = fn'. 

In order to relate the stress to other quantities through a matrix we wish 
to convert it into a single column matrix. We put f n — Xi , fa = X 2 , 
fzz = -X3 > /23 = /32 = -^4 , /31 = /13 = -^5 and /12 = /21 = X& . 

Changing to the X representation we find 

x' = aX' (8.1) 

where a is the matrix eq. (7.2). 
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Strain Theory 

If the dimensions of a body 'change, a point p = I p 2 J is moved to p -f- a v 

where tr 9 = I <rj J. A neighboring point p + •« is moved by an amount 
o-jh-u given by <r^„ = (Vo- c ) c m + <r P . The movement of p + « relative to 

/> is <T = (Tp+u — <T p = (VCc)cW- 

The 9 components of ( Vcr c ) c describe the sort of movement in the neighbor- 
hood of a point; they are the strain coefficients. If the strain matrix is 
e = (V<T c )c , a transformation x' = ax causes this to become aea e = {a\7a c a^ c 
and if aV = V and acr = </ so that a c a = <r c we have e' = ( V'<r e ) c if 

e' = oea c (8.2) 

When we arrange e as a single column matrix e we shall, following custom, 

take e 4 = — + — , e 6 = etc. This has the effect of moving the 2's of the a 
dxz dxa 

matrix to the conjugate position so that, while x transforms as x' = ax, 

e transforms as e' = aj e. 

We shall take tensions as positive stress elements, and elongations as 

positive strain elements. The shear strain, e c = (0, 0, 0, 0, e 6 ) becomes 

upon rotating through 45° about ;v 3 , e' c = ( ^ , — ^ , 0, 0, 0, J . This shows 

that to be consistent, a positive shear strain about x 3 must mean an expan- 
sion along the line Xi = x 2 and an equal contraction along the line Xi = —X2. 

A positive shear stress is one that tends to produce a positive shear 
strain. 

By superposing such strain elements we see that the e matrix (useful in 
displacement problems) may be formed from the e matrix (which is useful in 
stress strain relation) as 

(«1 , 2^6 2«b\ 
ie 6 e 2 \eA (8.3) 
\e h \e A e s / 

This slightly awkward relation is used solely to make the "work done in 
straining" expressible as 

2W = X c e = e c X (8.4) 

If the e's were taken as equal to the e's the work would be: 2W = Xtfi + 
X 2 e 2 + X 3 e 3 + 2Xie t + 2X b e b + 2X 6 e 6 . This would be awkward in some 
later problems. 
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If the scalar W is to be unaffected by a transformation a we must have 
W = e c X unaffected. If we write 

W = e<oT x aX = e^X' 

we have 

W' = W = e c X' 

if 

' -l 

e c = e c a 

when 

i = a?* (8.5) 

This substantiates our previous statement. 

Relation Between Stress and Strain 
If the strain in an elastic body is proportioned to the stress we may write: 

e\ = SuXi + S12X2 + ' " ' SieXe 
£2 = S21X1 + • • • 

Where the 5"s are elastic modulii. In matrix notation: 

e = SX (8.6) 

Conversely X = S~*e or if 5 _1 = C 

X = Ce (8.7) 

The C's are called elastic constants to distinguish them from the modulii S. 
As e = SX, a7 l e = aj l Sa~ l aX, and since a~ x e = e', (the representation 
of e on a new axis system related to the old one through the matrix a) and 
aX is X', then we may write (a7 1 e) = («7 SoT )(aX) as: 

e' = S'X' where S' = d?Sa* (8.8) 

Similarly operating on X = Ce we find 

X' = Ce' where C = aCa c (8.9) 

The energy required to cause the strain e is 

W - J X r de r = 1 X r e r = i Sr,X r X t (8.10) 

whence, if IF is a perfect differential, 

Srs m d *W = d * W =Str (8.U) 

dX r dXg dX g dX r 
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Similarly 



L" TS L/ a 



(8.12) 



This reduces the constants and modulii to 21 of each. 

If a transformation is performed that is permitted by the symmetry of 
the medium the elastic modulus matrix is unaltered. The monoclinic 
system has a binary axis. If we choose this as x 3 and rotate the axes 180° 

/-l 0\ 
about this by means of the matrix a = I — 1 I we have S' = 

\0 1/ 
a. Sac = S. 



100 










6116126 13O 146 1 BO 16 




100 





00 


010 


00 




61262262362462 fSu 




010 


00 


001 


00 




6136236336346 3 5636 




001 


00 


000 - 


-1 00 




140 240 34O 440 4 50 46 




000 -1 


00 


000 


0-10 




615625635645655656 




000 


-1 


000 


1, 




_ 160260 360 460 B6O 66 _ 




J) 


1 






Oil 012 °13 "J 14 Ol6 oi6 










Ol2 022 *->23 — »->24 «->25 626 










OJZ 623 633 634 — 635 0*36 










— Ol4 — 024 — 6 34 O44 645 — 'J 46 










— 6tf — 625 —S35 645 655 —Sm 










, *S"l6 




J 26 036 — 046 


Sti 


0*66 







Equating terms, those whose signs differ in 6 and 6' must vanish. 

Proceeding in this way through the 32 crystal classes we arrive at the 
ten following matrices that cover the elastic behaviour of all 32 classes. 



Triclinic System 
21 modulii 



6 = 



Monoclinic System 
#3 axis binary 
13 modulii 



6 = 



On 6)2 Ol3 OH Ol5 Ol6 

0]2 022 023 024 02 5 026 

•5*13 623 633 Su 636 636 

Su 6*24 Su Sa 645 Sit 

Sib S& 6*46 Sa Sbb 0^66 

Sie 626 Sw 646 0*66 6e6 

Sn Sn S13 6ie 

O]o 022 023 U 026 

o"i2 6*23 6*33 036 

S 4 4 6*45 

6*45 6* £5 

6*16 626 636 666 



The C matrix is 
entirely analogous 

(8.13) 



The C matrix is 
entirely analogous 

(8.14) 
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Phombis System 
x 3 binary 
9 modulii 



Tetragonal System 
x 3 a fourfold axis 
(Classes 9, 10, 13) 
7 modulii 



5 = 



5,i 5, 2 5,3 

5,2 5 22 5 23 

5 23 5 M 5 33 

S44 

5 M 

5 6 



The C matrix is 
entirely analogous 

(8-15) 



5 = 



5n 

5i2 

5 13 





5i6 



5.2 5i 3 
5n 5, 3 

5. 3 533 
5« 

5 4 4 
-5i« 



5,6 
— 5l6 








The C matrix is 
entirely analogous 

(8.16) 



Tetragonal System 
x 3 a fourfold axis 
Xi a twofold axis 
(Classes 11, 12, 
14, 15) 

6 modulii 

Trigonal System 
x 3 trigonal axis 
(Classes 16, 17) 

7 modulii ~ _ 



5 = 



5n 5, 2 5,3 
5iz 5n 5i 3 
5,3 5 13 5 33 
5 44 
544 
5 fl 



The C matrix is 
entirely analogous 

(8-17) 



5,i 5,2 5i 3 5i4 5,6 

5.2 5,1 5,3 — 5l4 5 2 5 

5.3 5,3 533 
0,4 — Oh O44 

■5 26 5 26 





25 26 
544 25i4 







25 2B 25 M 2(5,1 -5, 2 ) 



Trigonal System 
Xi trigonal axis 
:Vi binary (Classes 
18, 20, 21) 5 

6 modulii 
(alpha quartz) 

Hexagonal System 
^•3 a sixfold axis 
Xi a twofold axis 
(Classes 19, 22, 
23, 24, 25, 26, 27) 
5 modulii 



5n 5i2 5i 3 5m 

5,2 5n 5, 3 — 5w 
5 ]3 5i 3 5 33 
5,4 -5,4 544 
5 44 25m 
25 m 2(5n 



-5 J2 ) 



The C matrix 
is analogous 
except that 
C46 = C25 

Cb6 = C14 

Cee = 2 
(C„ -C12) 
(8-18) 

The C matrix 
is analogous 
except that 
Cm = Cu 
C K = h 
(C,i -C n ) 
(8.19) 



5 = 



5ii5, 2 5i 3 .0 
5i2 5n 5,3 
5i 3 5,3 533 _ a 
5 44 6-'. - 

0000 5 44 o:.--- 

2(5,1 -5,2) J 



The C matrix 
is analogous 
except that 

7 (Oii-C, 2 ) 
(8.20) 
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Cubic System 
Xi, X2 and x a 
fourfold axes 
3 modulii 



S = 



The C matrix is 
entirely analogous 

(8.21) 



Isotropic bodies 
2 modulii 

S = 



The C matrix is 
analogous except 
that 

^2 = 2 (Cll — Cu) 

(8.22) 



S u Sia S12O 
S n Su S12O 
Su Sn S n 
S u 
Su 
$44 

S11 5i2 5i2 ' 
Sn S u Su 

5i2 5l2 Su 

S 2 

s 2 
Jd s 2 

1S2 = 2 (5u —S12) 

Several Elastic Ratios in common use are given here for reference: 

Young's Modulus: A tension stress X divided by the component of 

X- 
strain in the direction of X , F,- = — - . If the coordinate axes are chosen so 

e» 

that the stress lies along X x , Yi = — . To find the value of Y in an arbi- 

011 

trary direction, (8, <p) find S' for a transformation that puts X' in the di- 
rection (0, <p) 

S' = aj Sot 
Where a is taken as form (21.4). Whence we obtain: 

f =— J = ClSaSu + s\siS22 + <htSu + S1S2C2SU + CiC25 2 555 + Ci5i52566 

+ 2C1S1C2S2SM + 2ciSiC2S2Si6 + 2C1S1C2S2SM ~\- 2C1S1S2S26 
+ 2CiSiS2Su + 2CiS2CiSt6 + 2CiSlS2C2Su + 2CiC2S2Su 
-\- 2C1S1S2S12 + 2SiC 2 S 2 .S23 + 25i52C2534 + 2ciSiC 2 S2Si6 

+ 2C1C2S2535 + 2C1S1C2S2S2& + 2S1C2S2SU (8.23) 

Rigidity Modulus: The shearing stress divided by the component of shear 
about the axis of shearing stress. For shear about Xi , 



Nx = 



1 



(8-24) 



Its value in another directions can be found as Y eip was above. 
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The bulk modulus: The change in volume per unit volume for unit hydro- 
static pressure is the bulk modulus, H. For a stress X c = (1, 1, 1, 0, 0, 0) 

e = (Sn -f- Sn + Sn , Sn + Sz3 + Saa , S31 + S23 + S33 , • ■ ■ ) 

£T = (d+ e 2 + € 3 = 5u+ S 22 + S 3 3+ 2^2+ 2S 31 + 2S 23 ) (8.25) 

This is obviously independent of the choice of axes. 

The Temperature Coefficient of the Elastic Modulii and Constants 
If 

C = C° + th + hi + hi + • • • (8.26)* 

and 

5 » 5° + tt + t*H l + / 3 i7 2 + • • • (8.27) 

(C° and S denote the values of the C's and S's for some standard tempera- 
ture / = 0) then as the transformations are 

C = a C a c and S' = aj SoT or 

C = a (C° + th + *V + t 3 h v • • -K 

and 

5' = 0^(5° + <ff+ Z 2 ^ 1 + tH 1 • • •) a -1 

we see that 

C = C 01 -f- «*' + W ■■■ (8.28) 

S f = S° l + /#' + /V • • • (8.29) 

where 

hi = a/?a c etc (8.30) 

H' = ct?Ha~ l etc (8.31) 

That is, the h's transform as the C's do, and the H's transform as the S's 
do. Consequently we may copy their respective forms from the C and S 
matrices for any particular crystal class. 

When the temperature coefficients of the constants or modulii are known 
in the form: 

Ca = C?/ (1 + tT Cii ) (8.32) 

Sn = Si, (1 + tT Bij ) (8.33) 

* The n of l n denotes the nth power of the scalar /; the n of /»" is merely another matrix, 
it does not mean a power. 
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we may write: 

h ti = C° i3 T Cij .(8.34) 

ffiy-^r^ (8.35) 

Multiplying (5) by (6) we get: 
SC = I = S°C° = S°C° + t(S°h + HC°) + t\S°ti + Hh + VC ) + ■ ■ ■ 

whence, for this relation to hold for all values of t: 

S°h + HC° = (8.36) 

whence 

h = -C°HC° 

(8.37) 

E = -S°hS° 

also 

S°l? + Hh + tfC = (8.38) 

so that 

h x = hS°h - C°H X C° 

(8.39) 

H 1 = HC°H - 5°A 1 5° 

From these we can compute the h's given the H's and vice versa. 

SECTION 9 

Temperature Expansion 

The change in the dimensions of a crystal caused by a temperature change 
can be considered as a strain. The shift of the terminus of a vector / 
relative to its origin is given from the strain matrix c by the equation Al = 
tel 

Since e is symmetric a proper choice of axes makes it possible to make the 
strain per degree a diagonal matrix, 

A4i o 0\ 

Al = tAl where A = [0 A 2 1 (9.1) 

\0 Aj 

As / and Al both transform as vectors, a transformation a causes A to 
transform as 

A' = aAcT 1 (9.2) 
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The elongation per unit length per degree in the direction (0, $) is 
/cos sin <p, sin 8 sin <p, cos tp 



Al Ba = 



whence 




cos v sin <p • • 
sin 8 sin<p • • 
cos = v 5 • \ 



Ale? = -<4i cos 2 sin 2 <p + A 2 sin 2 sin 2 <p-\- A 3 cos <p. 



(9.3) 



The strain can easily be extended to a function of t and / as follows: 



Al = tAl + ?BL 



(9.4) 



Applying the prefactor a to both sides and putting the idemfactor in 
between A and / and between B and / in the form / = a - a we have: 



aAl = t(aAa x )al + t 2 (aBa l )al or 
AV = t A'V + t 2 B'V where 
,4/ = aAaT 1 B' = o5o _1 



(9.5) 



SECTION 10 



Temperature Variation of the Isothermal Elastic Module: and 
Stress Variation of the Temperature Expansion Coefficients 

We can write the isothermal elastic modulus matrix at temperature 
8 + /as 

S* = S i0 + tH (10.1) 

and the coefficient' of temperature expansion at constant stress X as 

A = A° + LX .(10.2) 

Let us take a unit cube of crystal about the cycle indicated in the table; 
starting with the cube in the unstressed unstrained state at absolute tem- 
perature 8: 



Operation 


Change 
in Stress 


Change in Strain 


Temp. 


Heat at zero stress 

Apply X isothermally 

Cool at Const. X 




X 



-X 


tl° 

(5'° + M)X 

-t(A° + LX) 

-S i0 X 


9 to e + t 

6 + t 
+ ltoO 

e 







If we sum the strain changes in this cycle to zero we have 

H = L 
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so that we may write 

S = S + tH 

A = A° + HX 



(10.1) 
(10.3) 



This tells us that we may determine the temperature coefficients of the 
elastic modulii by measuring the effect of stress on the temperature ex- 
pansion coefficients. 

In a similar way we find that if the isothermal elastic constant matrix at 
temp. 6 + I is: 

C = C i0 + th (10.4) 

then the relation between temperature and stress at constant strain e is 

X = W (10.5) 

where 

B = B° + he (10.6) 

The Difference between the Specific Heats at Constant Stress and Constant 
Strain 

Writing for the specific heats at constant stress and at constant strain 
a v and o , respectively, we can perform the following cycle: 



Operation 


Change 
in Stress 


Change 

in Strain 


Temperature 


Work In 


Heat Out 


Entropy 


Heat at zero stress. . 
Restore zero strain 

Cool at zero strain . . 


- 


I A" 

-tA° 



e to e + 1 

e + t 
e + the 




p'.CA- 




—ptaP 

Q 
pt<r v 


—ptffP 

e + t/2 

Q 

e + t 

pta v 

•+i 



Equating the sum of the entropy changes to zero: 
Equating the work in to the heat out: 



( ff * _ a *) =°A° C CA° 
P 



(10.7) 
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Isothermal and Adiabatic Elastic Modulii 

Let us take a unit crystal cube at temperature 6, apply any stress X 
adiabatically, heat it to bring the temperature back to at constant stress 
then release it isothermally. The cycle is analyzed in the table: 



Operation 


U u 

Sat 

J5-C/J 

CJ 

X 


-X 


Change in Strain 


Temperature 


Work In 


Heat Out 


Entropy 
Change 


Apply X 
adiabati- 
cally 

Heat to 8 at 
const. X . 

Remove X 
isother- 
mally 


S"X 
t(A° + HX) 

-S i0 X 


6 to e - I 

e - ttoo 

9 


hx e s a x 
tX e (A° + HX) 

-±X c S io X 



— tpov 

Q 




—Ipffp 
9-t/2 

Q/9 



Summing the strains to zero: 

OS"' - S°)X = t(A° + HX) 

If we equate the total entropy change to zero we obtain an expression 
for Q that can be substituted in the relation "work in = Heat out." This 
gives us: 

-hX c (S io - S a )X + tX(A° + HX) - i ^ 
and from these two expressions we derive, writing <f> for S'° —S*: 



<t> = — p {A + HX)(A° + HX) C 

pa p 
which is, to the first order of the small quantities X: 





$ = JL A°A° C + 2HXA\ 

pa p 



and since X = Ce we have also 



-^= (A + 2HCe)A l 
paF 



(10.8) 



(10.9) 



(10.10) 



Whence we see that as the stress approaches zero as a limit <p approaches 

a 

0° = A° A° . If we write similarly C'° — C a = ^ we have multiplying 

p<r p 

S %0 = S° -f- <t> by C*° = C -f- ^ and dropping higher orders of small quan- 
tities: 

+ = -C io <pC io (10.11) 
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For example we have for quartz at 20° Centigrade 



293 X10" 



2.65 X 7.37 X 10" 



14.4 
14.4 

7:8 









(14.4 14.4 7.8 -0 0) 



3.13 3.13 1.69 

3.13 3.13 1.69 

1.69 1.69 .907 









10" 



and as 



S xo = 



10" 16 X 



1298 -166 -152 -431 



-166 1298 -152 

-152 -152 990 

-431 431 







431 


2005 
2005 862 
802 2928 



whence 0n = Sg (1 - -00241) 
0i2 = Sil (1 - .0189) 
0i3 = S\l (1 - .0111) 
3 3 = Siz (1 - .000917) 
For Rochelle Salt we have: 

4> = 293 X 10" 12 
1.79 X 15.5 X10 6 



38.0 24.2 28.5 

24.2 15.4 18.1 

28.5 18.1 21.3 









59.9 


38.1 


44.8 


4 





I o J 



(59.9 38.1 44.8 0) 



X 10" 
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As S t0 = 



4690 -795 -2180 

-795 3205 1691 

-2180 1691 2815 

60£ 































3060 








8020 J 



X 10" 



So that 



0i2 = S\° 2 (1 - .0305) 
<t>u = S[° 3 (1 - .0103) 
023 = S23 (1 - .0107) 



0ii = Su (1 - -0080) 
022 = S22 (1 - -0050) 
033 = Si° 3 (1 - .0076) 

044 = Su 
066 — ^66 
066 = S&6 

At the temperature of maximum piezo activity the components of 
for Rochelle are smaller by about 3£%. 

SECTION 11 

The Piezoelectric Effect 

Some crystals develop an electric charge when subjected to mechanical 
stresses. As far as the effect is linear it may be expressed by: 



Z/i — di\Xi -\- U12X2 ' • • "16-A6 

A = d 2 iXi + • • • 

Z> 8 = daiXi + • • ■ ^36-^6 



(lU) 



or in matrix notation 



D = dX. 



(11.2) 



where the 18 constants da are called piezo-electric constants, and D is the 
electric induction. 

On rotating the axes by means of a transformation a, the vector D be- 
comes D' where D' = aD. The stress transforms as X' = aX whence D = 
dX becomes D' = ada~ x X' or D' = d'X' where: 



d' = add 



(11-3) 




1 








1 








1 











-1 











-1 








1 
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If transformations permitted by symmetry are performed, the d matrix 
is unchanged. Class 3 has a binary axis only, if we choose this as x 3 and 
perform the transformation 

/-l 0\ 

o = ( -1 0] we find: 
V 1/ 



>dn dw d\2 d\\ d\$ 
d = i - 1 1 ( doi dn ^23 da dzb 

dsi da 0*33 0*34 dif, 



(—fl*ii —dvt —diz du dit — d^\ 
— dn —d-a —d 2 3 da d 2i —duj 
0*31 d 3 2 dzz —du —das <W 

For this to be consistent with the original d matrix the terms with conflicting 
signs must vanish. 
Applying similar analyses to each of the 32 classes we arrive at the set of 

matrices: 

Class 1 (asymmetric) /du dn dn du d i6 d i6 \ 

No symmetry d — I dn ^22 d<a d u d 2 b dn I (11.01) 

\dzi dsi dzz du 0*35 d&/ 

Class 2 (triclinic pinacoidal), center of symmetry d = (11.02) 

Class 3 (monoclinic sphenoidal A) du d u \ (sucrose) 

x* is binary d = lo 4*4*0 \ (11.03) 

\dn dn 4 das/ 

Class 4 (monoclinic domatic) Mi dn 4a ° ° 4*\ 

x z plane is plane of d = l dn dn dn d 26 I (11.04) 

symmetry \0 d 3i d 36 / 

Class 5 (monoclinic prismatic) center of symmetry, d = (11.05) 

Class 6 (Orthorhombic /0 d u \ (Rochelle) 

bisphenoidal) d = lo d 2b I (11.06) 

x h xt, x* binaiy . \0 dj 

Class 7 (Orthorhombic 

Pyramidal) /0 d u 0\ 

.-3 binary, *i and x 2 d = I du I (11.07) 
planes of symmetry Vsi ^32 </» 0/ 
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Class 8 (Orthorhombic bipyramidal), center of symmetry, d = 

Class 9 (Tetragonal 

bisphenoidal) d = 

x% is quaternary alternating 

Class 10 (Tetragonal 

pyramidal) d = 

x 3 is quaternary 

Class 11 (Tetragonal /0 

scalenohedral) d = I 

x 3 quaternary, xi and jcj \0 
binary 



(0 du d u N 

-</ 15 du 

Vrf 3 i -rfsi d*, 

/0 rf 14 Jib ' 

d xi -du 

\dn d 31 danO 0, 




d u 

du 

</36/ 

du 0^ 
-<* M 
0> 



Class 12 (Tetragonal 

trapezohedral) d = 

x 3 quaternary, xi and .r 2 
binary 

Class 13 (Tetragonal bipyramidal) center of symmetry, d = 

/0 d u 0\ 
du I 
V</ 3 i 4 4 0/ 



J = 



Class 14 (Ditetragonal 
pyramidal) 
x 3 quaternary 
xi and x% planes of symmetry 

Class 15 (Ditetragonal bipyramidal) center of symmetry, d = 

(dn —dn du du —Id-m 
— d& d-& Jig — du —2dn 
d 3l d zl 4 



Class 18 (Trigonal) 

trapezohedral) 
Xi trigonal, X\ binary 



d = 



Class 19 (Trigonal bipyramidal) 
Xi trigonal, x 3 plane of d = 

symmetry 

Class 20 (Ditrigonal pyramidal) 
x 3 trigonal, .r 2 plane of d — 

symmetry 



(11.08) 
(11.09) 

(11.10) 

(11.11) 

(11.12) 

(11.13) 
(11.14) 



Class 16 (Trinonal 
pyramidal) 
x 3 trigonal 

Class 17 (Trigonal rhombohedral) center of symmetry, d = (11.17) 




(d n -d n d u \ (Quartz) 

-du -2dn) (11.18) 

^0 / 

/ dn -dn -24\ 

-d 22 400 -24] (11.19) 

\ 0000 0/ 

/ du -2dn\ (tourma- 

-4 4 du ] line) 

\ 4 4 40 / (11.20) 



Class 21 (Ditrigonal scalenohedral) center of symmetry, d = (11.21) 
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Class 22 (Ditrigonal 

bipyramidal) 
x 3 trigonal, a* plane of d = ( -2d n ) (11.22) 

symmetry 
X2 plane of symmetry 

Class 23 (Hexagonal 

pyramidal) d = ( d u -d u J (11.23) 

x 3 Hexagonal 

Class 24 (Hexagonal 

trapezohedral) d = ( -4 0| (11.24) 

xa hexagonal, xi binary 

Class 25 (Hexagonal bipyramidal) center of symmetry, d = (11.25) 

Class 26 (Dihexagonal /0 4u 0\ 

pyramidal) d = lo <fe 6 OJ (11.26) 

x 3 hexagonal, x 2 plane \du d 3i 4 0/ 

Class 27 (Dihexagonal bipyramidal) center of symmetry, d = (11.27) 

Class 28 (Cubic tetrahedral- /0 <* M \ 

pentagonal-dedoca- rf=[0 <f M 1 (11.28) 
hedral) \0 dj 
x\ , X2 , x 3 binary 

Class 29 (Cubic pentagonal-icositetrahedral) d = (11.29) 

Class 30 (Cubic, dyakisdodecahedral) center of symmetry, d = (11 .30) 

Class 31 (Cubic, hexakis- /0 d u \ 

tetrahedral) <* = (o du 1 (11.31) 

Xi, Xz, x 3 quaternary \0 du/ 

alternating 
Class 32 (Cubic, hexakis-octahedral) center of symmetry, d = (11 .32) 

Whenever a center of symmetry exists the piezo-electric property vanishes 
since a center of symmetry requires d' = {—I)dl= —d= —d'. Also 
d = for the pentagonal icositetrahedral class. 

Classes 6, 11, 12, 24, 28 and 31 polarize only by shear. 

Classes 1, 3, 4, 7, 10, 14, 16, 20, 23, 26 can be polarized by hydrostatic 
pressure. As an example of this let us consider tourmaline (which is ditri- 
gonal pyramidal). For hydrostatic pressure, Xi = X 2 = X 3 , X 4 = X 6 
= Xq = 0, whence from the polarization stress matrices we find, A = 0, 
A = 0, A = (2d 3 i + d-x) X pressure. As d 3l - 0.75 X 10 -8 and d» — 5.8 
X 10~ 8 for tourmaline, we get 7.3 abcoulomlos per cm 2 per dyne per cm . 
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SECTION 12 

The Converse Peezo-electric Effect 

A stress X causes an electric induction 

D= dX (11.2) 

and a strain 

e = SX (8.7) 

If the charge is allowed to leak away a further strain occurs, at constant 
stress. This is the strain that would be gotten if the stress were originally 
applied with surfaces rendered conducting: 

e° = S°X (12.1) 

In the first sort of stress, the work per unit volume done on the crystal by 
establishing the stress X is: 

W = \X c e = \X C SX (8.4) 

The energy stored electrically in the medium is: 

W B = 2rDJt*D (12.2) 

while the work done on a conducting crystal is: 

W° = %X e S°X (12.3) 

If a crystal be stressed in its insulated state by expenditure of energy W, 
the charges then absorbed by an external circuit taking up energy We , 
the strain changes from e to e° at constant stress so that the stresses perform 
additional work 

W a = X c (e° - e) - X C (S° - S)X 

and the crystal is left containing energy W°. Whence 

W° = W - W B + W a (12.4) 

or: 

±X C S°X = \X C SX - 2ttD c 1T 1 D + X C (S° - S)X 

so that: 

X C (S° - S)X = lirDck^D 

If we substitute D = dX we find 

X C (S° - S)X = ArXcdck^dX 
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so that: 

S° - S = Airdjr x d (12.5) 

The change in strain caused by rendering the surfaces conducting is: 

e° - e= (S° - S)X = Awdck-'dX (12.6) 

If the crystal be now insulated and the stress removed, an induction of 
opposite sign will occur and because of the assumed linear dependence of D 
on X the new induction will be equal to the negative of the previous one. 
The induction D = —dX indicates an electric field: 

E = brk^D = brk^dX (12.7) 

Also, the strain will alter by an amount -e", where, since the action takes 
place with non-conducting surfaces: 

e" = SX 

This leaves a strain on the crystal, of amount: 

e' = e° - e" = (S° - S)X (12.8) 

From (12.6), (12.7) and (12.8) it follows that: 

e' = d c E (12.9) 

As the medium is in just the condition that an electric field E would put the 
unstressed medium, (12.9) is the equation of the converse piezo-electric 
effect. It is to be noted that the set of constants that relates polarization 
and stress is the conjugate of the set that relates electric field and strain. 
For convenience in notation the converse effect will be written as 

e= gE (12.10) 

where 

g = d c (12.11) 

Rewriting (13) as a~ l e = (a7 1 ga~ 1 )oE we see that 

e' = g'E' 

where 

g' = a?ga (12.12) 
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SECTION 13 
The Converse Piezoelectric Effect as a Non-Linear Function 

If the strain of a crystal is not strictly a linear function of the electric field 
causing it we must relate the components of strain to field terms of the sec- 
ond power as well as to first power terms. That is, the equation e = gE 
(which gives the strain e in terms of the electric field E through the 18 
constants g) must be modified to include terms £,•£, . All such terms are 
included in the symmetric matrix (E E c ). 

A transformation a that replaces E by a c E' also replaces E c by E c a so 
that (E E c ) is replaced by (a c E E c a), that is (E E c ) being self-conjugate, 
transforms similarly to the stress matrix. We may rearrange this as a one 
column matrix similar to the stress matrix X, as follows: 



(EE e ) = 



El 
El 
El 



E*E 3 

E&Ei 

{E&) 

We may now relate the strain to E and E through the two matrices g 
and G: 



= E 



(13.1) 



e = gE + GE 



(13.2) 



If transformations permitted by the symmetry of the crystal are per- 
formed, g' must equal g and G' must equal G, this allows us to simplify the 
matrices; g is no different than before and hence vanishes for all types having 
centers of symmetry (and for the pentagonal icositetrahedral class). 

Rewriting (1) as aj e = (aj ga~ )aE + etc Ga" aE we see that 



where 



e ' = g 'E' + G'E 



g' = a^ga 
G' = aT Ga 



(13.3) 



The matrix G transforms as the elastic modulii matrix does but G,y 5^ 
Ga . Applying G' = ajGa we arrive at the set of matrices that follow 
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Triclinic (36 consts) 

Gn Gn Gis Gu Gib Gie 
G21 G22 G23 G24 G26 G26 
G31 G32 G33 G34 G3B G36 
G41 G42 G43 G44 G45 G46 

G51 G62 G&3 G B 4 Gf,h Gb6 

Gei G62 Gg3 Ge4 G55 Geo, 
Orthorhombic (12 consts) 

Gn Gi 2 G, 3 
G21 Gn G 23 
G31 G 32 Gas 
G44 
G 86 
G 6 

Tetragonal Classes 11, 
12, 14, 15) (7 consts) 
Gn G12 G13 
d 2 Gn Gis 
G 31 G31 Gm 
G44 
G44 
G 6 



(13.4) 



(13.6) 



Monoclinic (20 consts) 
Gn G12 G13 G13 
G21 G22 G23 G26 
G31 G32 G33 G36 
G44 G 4B 
G M G 6B 
Gei Ge2 G63 Gm 

Tetragonal Classes 9, 10, 13) 
(10 consts) 

Gu G12 G13 G16 
Gn G, 3 -G n 
G31 G33 
G44 G 4B 
-G45G44 



(13.5) 



G12 
G31 




Gei-Gei 



(13.7) 



(13.8) 



Trigonal (Classes 16, 17) 
(10 consts) 

Gn G12G13 G14— G25 
G12 G11G13— G14 G26 
G n G31G33 

G41 — G41 G44 G4 B 2G B 2 
-G62 G52 — G45 G442G41 

2G 2B 2Gi 4 2(Gii-G 2 2)J 



(13.9) 





Trigonal (Classes 


18, 20, 21) 




(8 constants) 


Gu 


G12 G13 Gu 





G12 


Gn G13 — Gu 





G31 


G31 G 33 





G41 


-G41 G44 








G44 2G 4 i 





2G M 


2(Gn -Gn). 



(13.10) 



Trigonal (Classes 19, 22) 
(6 constants) 
Also Hexagonal (Classes 23, 24, 25, 26, 27; 

Gn G12 Gi 3 
G12 Gn G13 
G31 G31 G33 
G44 
G 44 
[0 2 (Gu - G U )J 



(13.11) 
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Cubic (Classes 28, 29, 30, 
31 and 32) (5 Constants) 

Gii G12 Gia 
G n Gii G12 
Gn G n Gn 
G44 
G44 
d 



(13.12) 



Isotropic Bodies 
(2 Constants) 
G11 Gi2 G12 ' 
Gia Gn Gn 
G12 G12 Gn 
G 
G 
G 
G = 2(Gn -Gn) 



(13.13) 



According to this analysis, all bodies suffer a change in dimensions when 
subjected to an electric field. These strains resulting from a field are 
generally much smaller than those strains e = gE present only in crystals 
lacking a center of symmetry. For example, quartz has a strain of about 
6.5 X 10~ cms/cm/ab volt. Glass in a field of 1000 practical volts per cm 
has a strain of about 4 X 10~ 12 , in a 100,000 volt field it has 4 X 10~ 8 . 
Rubber in the 1000 volt field strains by about 7 X 10~ 8 and in the 100,000 
volt field by about 7 X 10~ 4 . The 1st order quartz strain in these fields 
would be about 2.2 X 10~ 7 and 2.2 X 10~ 6 respectively. 

The Second Order Piezoelectric Effect 

If the induction stress relation is not strictly linear one can assume the 
induction to depend also on second order terms of the stress: 

D = dX + p (XX C ) 

where (XX C ) is a single column matrix formed from the 21 elements of 
XX C and p is a matrix of the 63 elements pn.i . . . p33,3- 

Since X transforms as X' = aX', (XX C ) transforms as X'X C = aXX c a c . 
In the same way that a was formed from a we can form a matrix a that 
transforms the single column matrix (XX C ) through (XX e )' = a(XX c Y- 

aD = adoT 1 X + ap(a)~ l a{XX c ) or 



where 



d' = d'x' + p'(xx c y 



d' = ada l and p' = ap(a)~ 



The first order effect is the same as before. With the relation p' = ap 
(a) -1 we could perform the operations of symmetry permitted by the 32 
crystal classes and obtain the reduced matrices. However since a has 484 
elements we shall limit ourselves to crystals with centers of symmetry. 

As X is unchanged by an inversion through the origin, a is the idemfactor 
for this transformation and a is —I, also (a) = /. Therefore D' = —D 
= D so that D vanishes. 
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It is to be noted that although there is a sort of reciprocity between the 
first order piezo effect and the converse effect, in that the matrices for one 
are the conjugates of the matrices of the other, there is no such reciprocity 
in the second order effects: if a center of symmetry exists no polarization 
can be brought about by stress either as a first order effect or as a second 
order effect; if a center of symmetry exists an electric field can cause a 
strain through the second order effect but not through the first order effect. 

Dielectric Constants at Constant Stress and at Constant Strain 

Let us consider a unit crystal cube, initially unstressed, unstrained and in 
zero electric field. We write k p and k v for the dielectric constant matrices 
at constant stress and constant strain, respectively, C B as the elastic con- 
stant matrix at constant electric field E, C° as the same for zero field. We 
study a cycle consisting of a strain caused by application of an electric field 
E at zero stress followed by a stress applied at constant E to reduce the 
strain to zero and completed by conducting away the electric charges at 
zero strain so that the body is left in its original state. The cycle is de- 
scribed by the table: 



Operation 


Change 
in Stress 


Change 
in Strain 


Change in 
Displacement Current 


Change 
in Field 


Energy Put In 


Apply E 

Apply — e. . . . 

Apply — E . . . 



-CH C E 

Cd c E 


d e E 
-d c E 




h k ' E 

4jt 


•J 

1 

4a- 


E c k"E 
Ecd C*d c E - 
E c (k" - k?)E 

E c krE 



whence C = C ; also 



k p - k v = 4tt dCd c 



(13.14) 



SECTION 14 

Pyro-Electricity 

If the electric polarization brought about by heating some kinds of crystals 
is simply a function of the uniform temperature change, that is if this 
polarization can be produced by taking the whole body quickly from the 
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uniform temperature 4 to the uniform temperature to + / the pyro-electric 
effect could be described by the equation: 



= t 



(14.1) 



where p is the pyro-electric matrix. 

This can be approached in another way by considering the polarization 
as due to the uniform strain. We may hence write, since X = Ce: (i.e . 
stress matrix = elastic constant matrix times the strain matrix) 

p = d X = dCe 

where e is the strain brought about by the temperature change /. If 

A4i \ 
^4 = I ^4 2 J is the temperature expansion matrix we have: 

\o o aJ 

[a; 

A, 

A, 




Now since d has 3 rows and the A matrix has but one column the product 
dCA has 3 rows and one column so that we may define p as dCA. 
As D of D = tp transforms by D' — aD, so does p: 

p' = ap 

When a center of symmetry exists a permitted transformation is a = —I, 
whence p = — p' = — p so that p = 0. No pyro-electric effect (on this 
theory) could exist for a crystal with a center of symmetry. 

If a binary axis exists and is chosen as x 3 we have 



P = tdC 




whence for this case 




(14.2) 



If another binary axis exists at right angles to this one we find p = 0. 
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This is seen seriously to limit the number of classes showing this kind of 
pyro-electric effect. In fact we find p = for classes 2, 5, 6, 8, 9, 11, 12, 
13, 15, 17, 18, 19, 21, 22, 24, 25, 27, 28, 29, 30, 31 and 32. The expression 

p = j ) describes the pyro-electric effect in the classes, 3, 7, 10, 14, 16, 20, 

h\ 

23 and 26, while the expression p = I p 2 I describes class 4, and only class 



v0 



1 is described by 




(14.3) 



It is to be noted then that this theory excludes many classes ordinarily 
described as pyro-electric, such crystals as quartz in Class 18 for example. 
Consequently it would seem that whether or not this effect exists we must 
seek elsewhere for the explanation of the effect in quartz. 

The effect can easily be explained as due to non-uniform temperature, 
which causes stress which in turn give rise to electric phenomena in piezo 
active crystals. For example a suddenly chilled crystal has its outer layers 
in a state of tension. This would produce just the pattern of positive and 
negative charges that one actually observes. As to whether the first effect 
exists, much argument between Lord Kelvin and others seems to have left 
the question still uncertain. 

In pyro-electric crystals we would expect to find a difference in the piezo 
constants measured isothermally or adiabatically. If a temperature 
change / causes an electric displacement D = pi the application of an electric 
field E should cause a temperature change t given by a relation such as: 

t = <pE (14.4) 

Also the temperature coefficient of expansion, A a (for a crystal with faces 
rendered conducting) would differ from the coefficient A { (for a crystal with 
an insulated surface). 

If a crystal at temperature h has suddenly applied to it a field E the 
temperature rises to to + <pE and the crystal strains, because of the converse 
piezo effect, by amount e = g a E where g a is the adiabatic converse piezo 
matrix. If the field is now removed isothermally a further strain gtE 
takes place. If the faces are short-circuited and the temperature restored 
to /o a further strain A s t = A 3 <pE takes place and the crystal is then in its 
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initial state. Equating the sum of the strains to zero we find (g a — gi) 
E = A ,<pE or 

g a - gi = A s <p (14.5) 

Let the initial state of a crystal be, temperature = t , stress, strain and 
field = 0. If the (electrically insulated) crystal is heated by amount /, 
a strain A it is caused and also an electric displacement D = pt. There now 
exists an electric field E = rrk~ pt. Let this field be discharged at con- 
stant temperature, giving a further strain of g,£ = ±-irgik~ l pt. The crystal 
is now short-circuited and if the initial temperature is restored a strain 
— A t t follows. The crystal is now in its initial state. If we equate the 
sum of the strains to zero we find : 

Ai- A = 4*giir l p (14.6) 

SECTION 15 

The Thermo-Electric Effect in Crystals 

It should be possible for an electric field to be set up by a temperature 
gradient. Let us assume that the vector T is the temperature gradient and 
is related to the vector field E through the matrix II by means of the equa- 
tion: 

/Tin n 12 n 13 \ 
E = TIT where II = I II 21 n 22 n 23 ) (15.1) 

\n 31 n 32 n M / 

Examination shows that II transforms through 

II' = aUa c (15.2) 

For Class 1 the II matrix has the 9 terms of (15.1). Class 2 has a center of 
symmetry. For a center of symmetry a = —I but a = —I causes no 
change in (15.2) so that class 2 has 9 constants. The thermo electro effect 
is not killed by the presence of a center of symmetry. The ordinary thermo- 
electric effect of metals is a case in point. 




If .v 3 is a binary axis a = I — 1 J and II reduces to 



/iin n 12 o \ 

n = fn 21 r^ o I (15.3) 

\o o uj 

Examination shows this form to answer for classes 3 and 4 and 5. 
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If xi and x 3 are binary (15.3) reduces to 

/an o o \ 

n = ( I^ ) (15.4) 

\o o iv 

which described classes 6, 7 and 8. 
For x 3 quarternary alternating: 

/ Hn nia o \ 

n = I -n^ n u o J (is.5) 

\ o o uj 

This is found to handle classes 9 and 10. If x 3 is quarternary alternating 
and Xi is binary: 

Aiu o o \ 

n= o n u o J (15.6) 

\o o iw 

which is found to cover classes 11, 12, 13, 14 and 15. For classes 16, 17, 
19, 23 and 25 II reduces to the form (15.5). 

If x 3 is trigonal and *i is binary the matrix is (15.5) which then handles 
cases 18, 20, 21, 22, 24, 26 and 27. 
For cubic crystals, not only are XiXz and .r 3 binary as for matrix (15.4) 
M /0 10\ 

but the vector 1 is an A 3 , for which a = I 1 I whence, for classes 

\l/ \l 0/ 

28, 29, 30 and 31 we find the matrix: 

Alii \ 

n = (o n,i o J (15.6) 

\o o n„/ 

Reports of a pyro electric effect in quartz should probably be attributed 
to nonuniform heating exciting the piezo electric effect. Reports of a pyro 
electric effect in such crystals as topaz and colemanite which have a center 
of symmetry and hence cannot be piezo electric should probably be at- 
tributed to this thermo electric effect. 

SECTION 16 
The Propagation of Light in Crystalline Media 

Maxwell's equations are: 

CV X B = 4*7 
CV X E = -B 



MATHEMATICS OF PHYSICAL PROPERTIES OF CRYSTALS 49 

when C is the velocity of light in free space, E is the (vector) electric field, 
j is the induction current and B is the magnetic induction. In a crystalline 
medium the current is given by Airj = k£ where k is the dielectric constant 
(matrix), whence: 

CV X B = k£ (16.1) 

CV XE= -B (16.2) 

As the divergence of the Curl is always zero: 

Vck£ = Vcj = 

(16.3) 

Vc-B = K J 

applying - to (16.1) and substituting (16.2) in the result: 

-C 2 V X V X E = kE or 

C 2 (VcV - VV,)£ = kE' (16.4) 

We shall try as a solution: 

V*""— ' (16.5) 



£ = *.Jt**cr-»t) 



where E is the vector amplitude of the electric field, i is y/—l, r is the 
radius vector from the origin to any point, q is a constant, n is the unit 
normal (at r) of surfaces of equal phase, and u is 2ir times the frequency of E. 
Substituting (16.5) in (16.4) we find: 

2 

E - n£ c n = ~- 2 kE (16.6) 

Examination of (16.5) shows that - is the phase velocity along n. Writ- 

2 



ing47r& j for £ and V 2 for — we have: 



? 2 

-1, ..,,-! V\ 



k j - nj e k n = —j (16.7) 



This equation is independent of the absolute value of j so let us restrict j 
to being a unit vector. 

Vck£ = = VckEoe il9 " cr - u,) = j c niqe i{qncT - ut) 

whence j e ti = (16.8) 

That is, the current is always normal to the direction of propagation. 
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Multiplying (16.7) thru by the prefactor /, and cancelling the term in 
j c n we have : 

i«*-V-5 < 16 - 9 ) 

This tells us that the velocity is a single valued function of the direction 

of the current. 

With the idemfactor I, (16.7) may be written: 

( -i v- V 1 

If we multiply this thru by \k - j 2 I J we get 

j = (k- 1 - ~ tV n{k- l j) c n (16.10) 

Multiplying this thru by n c and dropping the scalar factor (k~ j) c n: 

ne (k- 1 -^lj l n = (16.11) 

If the axes are so chosen that k is a diagonal matrix (16.9) and (16.11) 
become: 

T ,2 -2 .2 -2 

*£ + & + £ (16.12) 

w ' + __^ + *! — =0 (16.13) 

k n C* kn C 2 ^33 C 2 

Examination of (16.13) shows that (16.11) must have two values of 
V 2 for each value of the vector normal n. As V 2 is a single valued function 
of j there must be two distinct values of j (f and j" say) for any particular 
«• and given n, only waves having their current vectors in the directions 
of f and j" can be propagated. A ray in the direction N but not having 
its j in one of the directions/ or j" will be broken up into two components 
having their current vectors along;' and/' respectively. 

If the velocity Vi corresponds to j' and V 2 to j" we have by means of 
(16.10) since n' = «": 

j'j" = n c (iT 1 - ^ I)'' (k~ l - 2? ij 1 nik-'j'Uik-'Dcn 
(The quantities in the braces are scalar) 
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By means of the identity 

(m _1 - v" 1 ) = -u~\u - v)v~ x 
since 

the idemfactor can be multiplied into an adjacent matrix giving 

■•('-'-s<r('-'-9<r- 

=0-0=0 

so that/ and/' are mutually perpendicular. 

SECTION 17 

The Electro-Optic Effect 

The velocity of light in a crystalline medium is a single valued function 
of the unit current vector j 

7-i.^i (16-11) 

where c is the velocity of light in vacuo and k is the dielectric matrix, also 

dD 
j = D where D is — . 
di 

We developed the induction as a linear function of the electric field, 

deriving the relation: 

4ttZ> = kE (6.1) 

If the induction is not a linear function of the electric field we can improve 
on eq. (6.1) by adding second order terms: 

AvDi = k a Ei + k i2 E 2 + k i3 E 3 + k n El 4- kaEl + k*E\ 

+ \ huEvEz + \ fuiEsEi + . . . § hiiExEz 

or 

47rD,- = (ka + h a Ei + 1 h i9 Et + h haE*) Ei 

+ (ft« + § ft«& + h i2 E 2 + I *,«Eb) -^2 

+ (k a + i Aiftfi, + \ huEi + A^s) E 3 (17.1) 
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Examination of (17.1) suggests that we might consider the k's as being 
lineraly modified by the field. Writing Hsa single column matrix: 



K = 



feu 




\kA 


£22 




K 2 


£33 







hz 






k\Z 






[kn\ 




kJ 



we may write 



K = K° + hE. 



(17.2) 



where K" is the dielectric matrix for vanishingly small fields. 

We can develop the modified reciprocal matrix in the same manner: 



K~ l = 





ku 


ku 


few 1 ] 


k- l = 


k\2 


kn 


£23 




k\3 


few 


few 



(KT 1 ] 




ku 


Ki l 




&22 




= 


kit 

kiz 


[*?] 




feji 



where 



K~ l = K' 1 + zE 



(17.3) 



It is to be noted that K~ l is not the reciprocal of K but merely a symbol 
for the single column matrix formed from k~ l in the usual way. Taking 
reciprocals of both sides of k' = aka c we find (& -1 )' = ak~ a c . That is, 
JT 1 transforms exactly as did k. Whence, K~ l transforms exactly as K 
does, i. e. 

FT 1 ' = *K- 1 

We can rewrite (17.3) as 

aRT 1 = a K' 1 ' + (a za c ) aE 



or 



K7 V = {K~ i y + z'E' 



z = a za c 



(17.4) 



In case a l = ac the z's transform as do the conjugates of the piezo- 
electric constants, d. Of the transformations permitted by the symmetry 
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of the 32 crystal classes only those of the trigonal and hexagonal systems 
fail to have a~ = a c . These 12 classes must be examined individually but 
the other classes may have their z matrices copied from the corresponding 
d c matrices. 
Applying z' = aza c for a rotation of 120° about X* we find for class 16 



z = 



Zll 


— Z22 2]3 


— Zii 


Z22 Zl3 





Z33 


=41 


zu 


ZB1 


-Z41 


— Zn 


-zu J 



The remaining 11 classes may be derived from class 16 by operations for 
which either a -1 = a c or a center of symmetry exists. Consequently, we 
may form our z matrices from the d,'s in all cased if we leave out the 2's. 
The electro-optic effect can be put in terms of the polarization instead of 
the field by substituting in (17.3). 



E = 4tt (k-I)' 1 P whence 
KT l = K' 1 ' + V P 



(17.5) 



where 



Conversely 



V = 47rz(fc - I)' 1 



(17.6) 



The 77 matrices transform exactly as did the z's and hence may be formed 
from the d c 's but omitting the 2's. 

SECTION 18 
The Piezo-Optical Effect 



If the dielectric constants of a crystal are changed by the application of 
stress, this may be represented by: 

- 1 = (ir 1 ) *** 



where the 36 constants m • • • tm are stress-optical constants. 
We may then form k~* as 



k~ l = 



'kit + TTliXi, 

£12 + irt*Xi, 

\&13 4" TnXi t 



ku + 7T6iXf, 
^22 + TtiXii 

kn + TTaXi, 



k^ + T 6i X^ 
^23 + TTaXi 
^33 + TftiXiJ 



(18.1) 



(18.2) 
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As the velocity of a light ray of unit current vector j is given by 



Jc* 3 2 



(16.9) 



We can, by (18.2) and (16.9), compute the change in the velocity caused 
by the stress, if we know the constants t. 

Altering (18.1) to aK" 1 = aKT 1 ' + aira~ l aX we see that: 



K 



—v 



= K~ V + t'X' where *> = aicct l (18.3) 

The alteration of K' 1 can be expressed as a function of the strain by 
substituting ce for X in (18.1). 



KT 1 = K' 



m = ire 



+ iree = K x " + me 
t = tns 



(18.4) 
(18.5) 



Operating in (18.4) as we did on (18.1) we find m transforms as 

m' = ama c (18.6) 

Applying the crystal symmetry operation to these matrices shows that 
they reduce to the following 



Triclinic system 
36 constants 



Monoclinic system 
Xt is binary 
20 constants 



Orthorhombic 
system gg 
is binary 
12 constants 



Til 7T12 7T13 Tl4 Tl5 7Tl6 

7T21 7T22 ^"23 7T24 T>5 T26 

7T31 7T32 T33 T34 7T3B T36 

T41 T42 7T43 T44 7T45 ^"46 

T51 7T52 '""53 7TB4 T66 7TB6 

7T61 T62 7T63 T64 T"65 ^"66 

T11 7T12 T13 Tie 
T21 7T22 7T23 7T26 
T31 7T32 7T33 T36 
T44 7T45 
7T54 7T65 
7T61 7T62 7T63 T66 

Til T12 T13 

T21 T22 T23 

T31 T32 7T33 

T44 

T 6 6 

T 66 



The m matrix 
is entirely 
analogous 



(18.7) 



The m matrix 
is entirely 
analogous 

(18.8) 



The m matrix 
is entirely 
analogous 
(Rochelle salt) 
(18.9) 
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Tetragonal system 
x> is a four-fold 

axis 
(Classes 9, 10 

& 13) . 
9 Constants 

Tetragonal system 
%i is a four-fold axis 
Xi is a binary axis 
(Classes 11, 12, 14 & 

15) 
7 Constants 



7Tn T12 Xi3 7Ti6 

1Ti2 Tn 7T13 — 7Ti6 

7T31 7T31 T33 

7T44 

7T44 

7T61 — 7T61 7T66 



The m matrix 
is entirely 
analogous 



(18.10) 



ITu 7T12 7T13 

7Ti2 7Tn 7Ti3 

7T31 7T31 1T33 

7T44 

7T44 

7T66 



The m matrix is entirely 
analogous 



Trigonal system 
Xz is a trigonal 

axis 
(Classes 16 & 

17) 
11 Constants 



7Tn 7T12 7T]3 T14 — 7T25 

7T12 7Tn 7T13 — T14 7T26 

7T31 T31 7T33 

7T41 7T41 7T44 7T45 2^62 

TB2 7T52 — 7T46 7T44 2tTh 
7T25 7T14 (tu — T12) 



(18.11) 

The m matrix is analo- 
gous except that 
w 46 = m 62 

W 6 6 = W41 

mn — W12 
(18.12) 



Trigonal system 
#8 is a trigonal 

axis 
%i is a binary axis 
(Classes 18, 20 & 

21) 
8 Constants 

Hexagonal system 
Xz is a sixfold axis 
X\ is a binary axis 
(Classes 19, 22, 23, 
24, 25, 26 & 27) 
6 Constants 

Cubic system 
3 Constants 



7Tn 7Q2 Tl3 7Tl4 

7I"12 7Tn 7I"13 — 7Ti4 

7T31 7T31 7T33 

7T41 — T4l T44 

7T44 2X41 

7T14 (Tin 



7Tn 7T12 7Ti3 
7Ti2 7Tii 7T13 
7T31 7T31 7T33 
7T 44 
7T44 

[00000 (tu 



— 7Tl 2 ) 



The m matrix is analo- 
gous except that 
Tom = W41 

wn — m n 



ntw = 
(quartz) 



(18.13) 



— 71-12) 



The m matrix is ana- 
logous except that 

mn — mi2 

mm = 2— 

(18.14) 



7Tll 7Ti2 7T12 

7T12 7Tu 7Ti2 

7T12 7T12 Tu 

7T44 

7T44 

7T44 



The m matrix is en- 
tirely analogous 



(18.15) 
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For isotropic bodices, the t matrix is formed by setting x« = (*fc— *n) 
in the tt matrix of the cubic system; the m matrix is similarly formed by 

putting w 4 4 = ~ • 



Isotropic 
bodies 



TTn 7T12 7T12 

7Ti2 T11 ""12 

7Tl2 7Tl2 TTn 

irn-7ri2 

71-11-71-12 

7TH-T12J 



The m matrix has 
ma etc = - 

(18.16) 



SECTION 19 
Application of the Electro and Piezo Optical Effect , 
In the equations K~ l = IT 1 ' + zE and K' 1 = K~ l ° + mE, etc. the K~ 



are to be used in forming k 1 for the equation giving the velocity of the light 

used namely— = jck~ l j. Obviously then K 1 should be formed from the 

squares of the reciprocals of the refractive indices, the lower three members 
being zero. After applying the electric field or strain a transformation of 
coordinates may be necessary to rediagonalize, i.e. make KJ = KT = 
KJ 1 ' = 0. From the rediagonalized K~ l we may write the new principle 
refractive indices by taking the reciprocals of the square roots of KT , 



K2 1 ' and K z 



—V 



It should be noted that if K { = K { + A, then 

M<3 



n't + Am = n'i - ^ Ai 



(19.1) 



For a given direction of the wave normal there are two velocities, a wave 
splitting into two components traveling with different velocities. By defini- 
tion the refractive indices, /*« and n b in a given direction are the normal 
velocities in that direction divided into the velocity of light in free space. 

Ilia . , lltb 

Whence in a path length I there are — waves in one component and — 

waves in the other, where X is the wave-length in vacuum. Consequently if 
I is the thickness of the crystal along that path the two components can 
recombine after passing through the medium but they are out of phase by 

y = i ( Mo - fx b ) whole waves so that the light which entered as plane polar- 

A 

ized will leave elliptically or circularly polarized, except when 2t is an integer. 



MATHEMATICS OF PHYSICAL PROPERTIES OF CRYSTALS 57 

The quantity B = p a — /i& is known as the birefringence. 

7=^ (192) 

If a phase difference of -fa wave can be just detected, using a wave length of 
6000 A and a path length / = 1 cm the just detectable birefringence would 

be B = 1- = 2 X 10" 6 ; if the path were 10 cms the detectable B would be 

I 

0.2 X 10~ 6 . Obviously this detectable difference between refractive indices 
is much smaller than could be detected by measuring each refractive index 
and subtracting. 

It is customary to choose the coordinate system so that looking along x 2 
the very lowest refractive index is for polarization in the plane of xi and the 
very highest for polarization in the plane of x 3 . That is, the x 2 axis is the 
axis along which light should be passed to get the greatest birefringence. 

Birefringence in any Direction 

If the axes are so chosen that K is diagonal and K 3 > K 2 > K\ then, 
somewhere in the plane perpendicular to x% are two directions, the optic 
axes, along which there is a single normal velocity. These directions make* 
equal angles V with the xa axis where 



sin V = ^V KI 1 - KT l 

„ = (19.3) 

cos 



. /X-2 — Kz 



Also the two refractive indices n a and n b for a wave normal making angles gi 
and g 2 with these optic axes satisfy the equation: 

\ = (KT 1 + KT 1 ) + (KT 1 - KJ 1 ) cos ( gl - g 2 ) 

Ma 

\ = (KI 1 + K?) + (KI 1 - KI 1 ) cos ( gl + g % ) 
Mb 



whence 

2 






l\ 2(nb — Ma) (Ma + M&) 



2 / 2 2 

Mb/ MaM6 



. = (KI 1 - K s x ) cos (gi - g 2 ) - cos (gi + g 2 ) 
* Theory of Optics, P. Drude, pg. 320. 
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as B = nb — n a is the birefringence we have: 

B = -^p- (-KT 1 - KJ 1 ) sin gl sin g 2 

Ma + M6 



(19.4) 



(19.5) 



By spherical trigonometry: 

cos gi = cos V cos 6 + sin V sin 6 cos <f> 

cos g 2 = cos V cos 5 — sin V sin cos 

where is the angle the wave normal makes with x 3 and </> is the angle the 
plane containing the normal and x 3 makes with Xi. 

From (19.5) it follows that: 
sin gi sin g 2 = 



• 



(1 — cos 2 V cos 2 d — sin 2 V sin 2 6 cos 2 cf>) 2 



— 4 sin V cos F sin 6 cos" cos 



(19.6) 



Hence if the rediagonalized K is 

to 1 " + Ai 

K7 1 ' + A 2 

i^" 10 + A 3 





then 



5 = 



2 2 
fia fib 

Ha — lib 




(KT 1 ' - Kz" 1 ' + Ai - A 8 ) - {K?° - KJ 1 ' + A 2 - A 3 ) cos 2 6 
- (Kx l ° - Kz l ° + Ai - A 2 ) sin 2 d cos 2 <£ 2 - l(KT l " - KJ 1 ' 

+ Ai - A 2 )(£7 r - XT 1 ' + A 2 - A 3 ) sin 2 cos 2 cos 2 . . (19.7) 



5 = Bo + y (AiA 3 ) sin gi sin g 2 



(3') 



For most practical purposes we may take 

2 2 3 

fig fib _ M_ 

Ma + M6 2 

where /x is some intermediate value of the refractive index.* 

2 

* Note: It might seem that as K' 1 = K7 10 + A,- gives us ai,- = ti? — ^ A,- + fM?A? • • • 

we could form the 3 principal birefringences directly from the fi's instead of using (6a • • • e). 
From the n expressions we would get B3 = /x 2 — m — KmIA: — m?Ai) + f which differs 
from 6b if fit ,£ m . Equation 66 is correct; the one from the fi expression is an approxi- 
mation. 
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In a few special cases (18.7) may be simplified. If (0, 0) falls along an 
optic axis gi = V — V° and #> = V + V° whence 

B A = fL. [(Ax - A 2 ) cos 2 V° - (A 2 - A3) sin 2 V°] . . . . (19.61) 



if = 



if = 90 c 



J5 = ( M3 - mi)(1 - cos' sin 7°) 



2 2 

i? = M2 - Mi + -^- (Ai - A 2 ) 
Ml + M2 



2 T7-ON 



(19.62) 



+ |- (Ax - A 3 - (A x - A 2 ) cos 2 0) . . (19.63) 



if = 



B = (Ma - Mi) (sin 2 V° - sin 0) + £- (A x - A 2 - (Ax - A 3 ) sin 0) . . (19.64) 

if = 90° 

B = (jus - Mi)(l - cos 2 cos 2 F°) 

+ \ (Ax - A 3 - (A 2 - A 3 ) cos 0) . . (19.65) 

77*e Electro Optics of Quartz 
For quartz, in the equation X - = K~ + z£ 

211 ol 



-2x1 



Z - 2 41 

-z« 

-zix y 

Obviously the £3 component produces no effect so we shall examine the 

effects due to the components £x and E 2 separately. If E = I I Ei 

W 



If E is in practical volts* 
Z11 = 0.47 X 10~ 10 
Zu = 0.20 X 10" 10 



IT 1 = 



KT 


I + 2u£l ] 


KT 


° - znEi 




&° 




ZnEi 













which can be diagonalized 



* Computed from F. Pockels data, see his Lehrbuch der Kristall-Optik, (B. G. Tuebner). 
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60 



BELL SYSTEM TECHNICAL JOURNAL 



— Z41 Ei 



by a small rotation about x x of amount 8 = =1 ~i giving 

Ai — A3 



A" 1 - 



AT 1 


+ 


Zll-Ell 


AT 1 


— 


Z11E1 




AT 1 
























, 



3 
Mi - j X 0.47 X 10" 10 £i 

and n = mi + ~ X 0.47 X 10" 10 £, 

M3 



The greatest "added birefringence" is gotten by viewing along Xs , when 
AB = 1.544 3 X 0.47 X 10~ 10 Ex . If £1 = 10 4 A5 = 1.73 X 10" 6 a quantity 
detectable if the path length is about 1 cm. Viewing along x 3 (the optic 
axis) is complicated by the rotation of the plane of polarization in quartz. 
Homogeneous strains have never been found to alter this rotation, but the 
rotation complicates and partly masks the birefringence phenomena. If 
M 



E = 



£2 we find 



AT 1 " 

at 1 ; 

j XT 

A - Q 

— Z41.E2 
— Z11E2 

Rotating the coordinate axes through 45° about x 3 then applying the 
transformation 



a = 




where 



we find: 



i" 



KT l " = 



a = 



AT — A~3 



aT 1 
aT 1 


+ Z11E2 
— Z11E2 
A3 


















3 

Mi — -x-ZuE* 
3 
M = Ml + y z ii -^2 

M3 
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which is identical to the n for field along x t , but the final axes in this case 
do not coincide with the final axes for E = E\ , but again, the greatest added 
birefringence is utilized by viewing along x 3 . In the second case the Nicols 
would be best set along Xi and x 2 , i.e., at 45° to x[ and x 2 whereas in the 
first case they would be best set at 45° to xi and x 2 . 



The strain Optics of Quartz 
K~ l = IT 1 ' + me, 



m = 



Wn W12 1Yl\Z 


Wl4 


Wi2 Wn W13 


— W14 


W31 W31 W33 





W41 W41 


W44 





W44 







W14 









W41 
mn — W12 





tnu = 


.138 




nin = 


.250 




W13 = 


.259 


where* 


niu = 


.029 




nisi = 


.258 




ni33 = 


.098 




ma = 


-.042 




mu = 


-.0685 



If the strain is a simple tension along Xi , 



K' 1 = 



KT 1 ' 


+ f»uei| 


KT V 


+ Wl2Ci 


ST 1 ' 


+ W 4 l€i 




W 4 iCi 



J 



which diagonalizes, 
through a small 
transformation to: 



IT 1 ' = 



k; 1 ' 


+ 


*»n«il 


k; 1 ' 


+ 


Wi2 6i 


ej 10 


+ 




W 3 i€i 







t 



applying 18.63 or 18.64 we find the birefringence along X\ to be: 

3 
Bu = Ms - Mi + j (W12 - m»i) «i = .0091 - .0148 «i 

Similarly the birefringence along x 2 is 5 i2 = .0091 — .225 £1 and B\z = 
— .207 £1. With a strain of 10 -4 , which is about a tenth of the breaking 
strain, B IZ would be 20.7 X 10 -6 , a quantity detectable in a thickness of one 
millimeter. 

The values of B21 ••• B i3 corresponding to birefringence along tf 6 for a 
tension along x 2 etc., can be computed in just the same way. But B 6 i • ■ ■ 
^63 require rotations of 45° about x 3 to diagonalize, so the birefringences 
can be computed by setting 6 = 45° in equations (19.6). 



* Lehrbuch der Kristalloptik— F. Pockels. 
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The following table summarizes these simple strain birefringence effects, 
the rows indicating the strain and the columns the direction of light passage. 



Bo - .0148 ei 


Bo - .222 ei 


.207 ex 


Bo + .222 ft 


Bo - .0148 ez 


.207 eg 


B - .298 ez 


Bo - .298 ej 





Bo + .0536 ^ 


Bo - .0536 ei 


-.107 ej 


(5o+0 


Bo + 


-.107 «b) 


CSo + 


Bo + 


.208 e t ) 



(the parentheses indicate the 45° transformations). 
A similar table for the electro-optic effect in quartz is 



B - .87 X 10" 

Bo+0 

Bo + 



% 



Bo + .87 X 10' 

-Bo+0 

-Bo + 



°£i 



1.74 X 10" 10 £i 




Since a driving voltage of E\ = 100 volts may, due to the building up of 
oscillations, cause a periodic strain of e = 10 _ in a quartz plate, it would 
seem from the foregoing that 99.99% of any birefringence change must be 
due to the mechanical effect. 



The 18° Cut Crystal 

A crystal, the thickness of which is along the electric axis, xi , the width 
making an angle 6' = 18° with the optic axis, #3 , can be caused to oscillate 
with a simple motion along its length. (If 6' is not about 18° or 72° the 
oscillation is not a simple extension along the length, as is shown by the node 
which then lies diagonally across the crystal.) On a set of axes defined by 
the edges of the crystal block, x\ being in the direction of the thickness or 
Xi, x'z in the direction of the width and makes an angle 6' with # 3 , #2 is the 
length and makes an angle 6' with x 2 ; on these block axes the strain is £2. 
Rotating the axes about x[ through an angle 6' we find the strain expressed 
on the crystal axes to be: 



-1 / 
e = a e = 










N 


2 
COS 


d' 






■ 2 

sin 


e' 






sin 


6' 


COS 


e' 


















J 



ei 
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whence 

K~ l — K~ l + we gives us 

f R~[ l ° + (W12 cos 2 0' + W13 sin 2 5' — w i4 sin 0' cos 0')?? | 
Xr 1 " + (mu cos 2 0' + W13 sin 2 0' + « M sin 0' cos d')ej 
KJ 1 ' + (w 3 i cos 2 0' + fw 33 sin 2 B')ei 

— (w 4 i cos 2 0' + m»A sin 5' cos p')«2 







ir 1 = 



A small transformation removes the 4th term without altering the others to 
the first power of small quantities. 

To obtain the birefringence along the width x 3 , we set 8 = 8' in equation 
(18.65): 



Bv = {fi3 — Mi) sin 2 8' + -=■ \ mn cos 2 8' — Wu cos 4 0' + {m a — W33) sin 4 8' 

. 2 n , 1 + cos 2 6' . 2 „, 2 -A / 

— m u sin 0' — w 3 i sm 0' cos 0' > £2 

which, for 8' = 18° is B 3 = .00087 + .20 d 

For the birefringence along the length x 2 _ we set 

= 90° + 0' in (6«) giving: 

3 
B = (n3 - Mi) cos 2 0' + ~ {W12 cos 2 0' - m n cos 4 0' + (w^ — m u — w 33 ) 

sin 2 0' cos 2 0' - % sin 0' cos 0'(1 + sin 2 0') }d 
which, for 8' = 18° is £ 2 = .00824 + .049 ei 

SECTION 20 
Transverse Isotropy 

A material that has identical properties in all directions normal to a given 
line is called transversely isotropic. Any line parallel to this line may be 
considered as an axis of transverse isotropy. 

Dielectric Properties, Optical Properties, Thermal Expansion 

With respect to these, a transversely isotropic material behaves as does a 
uniaxial crystal, only two constants being needed to describe each. For 
example, the displacement current in terms of the electric field and the dielec- 
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trie constant matrix is D = — kE where, if x 3 is the axis of transverse 
isotropy: 



D = 




(20.1) 



Elasticity 

We must find the forms of S and C, (the elastic modulus and elastic con- 
stant matrices) that are not changed by rotations about the axis of trans- 
verse isotropy. We can simplify the work by starting with the crystal 
class that has hexagonal symmetry only. On applying the transformation 
S' = a c Sa = S for arbitrary rotations about X 3 we find no further simplifi- 
cation follows. Hence the 5 and C matrices can be copied from those for 
the Hexagonal Pyramidal Class. 

The Piezo-Electric Effect 

Again choosing x 3 as the axis of transverse isotropy and starting with 
hexagonal symmetry about X 3 we find that in order to be invariant to all 
rotations about X 3 the matrix must simplify to: 



/0 0> 

d = [0 000 

W d 31 d 33 0; 



(20.2) 



A pitch solidified in an electric field would probably exhibit this kind of 
piezo electric behaviour. It might also be expected to show an electro optic 
effect governed by a matrix like the conjugate of the above matrix. 

SECTION 21 
Appendix 

Transformations 

A counterclockwise rotation of the axes through an angle <f> about the *i 
axis is represented by the matrices a and a as follows (where c is written for 
cos <f> and s for sin <£) : 




a = 



1 





0] 







c 2 5 2 2sc 










s 2 c -2sc 










2 „ 2 r\ 
—sc sc c —s 










c 


—5 







5 


c 





(21.1) 
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A counterclockwise rotation about Xt is given by: 




c 2 


2 
5 


-2sc 


1 








s 2 


2 

c 


2sc 








C 5 


sc 


— sc 


c 2 -s 2 








-5 c 



(21.2) 



A counterclockwise rotation about x» is given by: 



a = 




2 2^ 
C S 








2cs 


s 2 c 2 








— 2cs 


1 














c 


—s 








—s 


c 






— cs cs c 2 — 5 2 



(21.3) 



In case one wants only the value of a tensor property in a given direction 
not all the elements of a and a need be used, but only a row or column. A 
special case is that of computing such a property in the direction (6, cf>) 
of polar coordinates. The X\ axis is chosen in this direction; Xt and x 3 are 
not determined. Writing C\ for cos 9, c- 2 for cos <f>, s l for sin 9 and s 2 for sin <f> 
the required matrices are 




2 2 
ClS2 

2 2 
SlS 2 

2 
Cl 

S1S2C2 

Ci Ci Si 
C\S\S% 



(21.4) 



From these the (11) term can be computed for any tensor. 

A few special transformations needed constantly are: A rotation of 180 c 
about #g: 




1 












1 












1 





-1 






-1 
, a = 


= a 








-1 












1 







(21.5: 
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A rotation of 90° about x 3 : 

10 
10 
10 
- i 
1 


A reflection in the plane perpendicular to #3: 







0] 














1 














-1 

4 



= a c ..(21.6) 



/10 





= 01 





\0 


-1 



1 

1 

1 













-1 






6) 











-1 


i y 



= a c ..(21.7) 



A cyclic interchange where .r 2 replaces x\ , etc. : 

(The line making equal angles with x h x 2 , and x 3 is a three-fold axis) 




1 ol 


10 


10 


1 


10 


10, 



(21.8) 



a = 



-1 
a = oc, 



(21.9) 



A cyclic interchange, where —x 3 replaces x h etc.: 

(The line making equal angles with x h .v 2 , and x 3 is a six-fold axis of the 

second sort) 

(0 1 o' 
10 
10 
a ~~ 10 
1 
1 0, 

A rotation of 90° about x 3 combined with a reflection in the x 3 plane: 



a~ 1 = a c .(21.10) 





1 








1 








1 











1 








-1 






0-1 
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A rotation of 60° about x 3 : 



1 V3 



2 

V3 

2 





1 



To form a ' substitute 
— \/3 for \/3 in a. 



o 1 









V3 V3 






1 

2 

V3 

2 










2 
1 

2 



V3] 
2 

V3 
2 








•-i 



.(21.11) 



A rotation of 60° about .r 3 combined with a reflection in x>: 



a = 



To form a 1 substitute 
— y/3 for -\/3 in a. 









1 

4 


!° 





1 y/l 

2 2 







3 
4 



1 






\/3 1 
2 2 





a = 








1 V3 

2 2 





-1 







- 


\/3 1 
2 , 2 



V3 V3 



V3 
2 

V3 
2 






-r 



..(21.12) 



A rotation of 120° about .r 3 : 



a = 



1 y/l 



2 

V3 



1 



To form a ' substitute 
— V 3 for \/3 in a. 



1 


7° 








4 


4 






3 


1 ^ 








- 








4 


4 









1 














1 


_V3 






2 


2 








V3 


_ 1 






2 


2 


V3 


-V|o 








4 


4 







V3l 
2 

V3 
2 







(21.13) 
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Inversion through the origin (a center of symmetry) 




a = 



A rotation of 180° about %i'. 




1 o ol 


10 


10 


10 


10 


1, 



= /. 



(21.14) 



-l 

a — a. 









0-1 



10 





10 





10 





1 








-1 



= a r ...(21.15) 



A reflection in *i plane: 




A reflection in the x 2 plane: 




10 
10 
10 
1 
0-1 




1 
1 
1 









0-10 
1 












-1 









0-1 



-1 

a = a, 



..(21.16) 



-1 
a = a = a c . 



(21.17) 



In computing the electro-optic and mechanico-optic effects we need a 
transformation that will restore to diagonality a matrix that has very small 
but symmetrical off diagonal terms. This transformation we call a small 
transformation. Such a transformation has its matrix differing but slightly 
from an idemfactor. 



If 



k = 




A12 A 3 i > 

&22 ~T" A02 A23 

A 2 3 #33 + A 3V 
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we assume that it can be rediagonalized to the matrix fe' by means of the 

/l Sis 8»\ 

transformation: fe' = 8kd c where S — I 621 1 fet J 

VSI 532 1 / 

Since 8 n S 2 i + 5i 2 5n + 5ia 5 2 a — we find that, to the first order of small 
quantities 5,7 = — 5>,. 

Expanding 8k~8 r to the first order of small quantities and equating the non- 
diagonal terms to zero we find that: 



5l2 — 



823 — 



A; 3 



£11 "~ kn £22 - £33 

Therefore, to the first order of small quantities A ,7-: 



d 5,, = 



£33 — £11 



(fen + An A12 A 8 i \ 

A12 fe 22 + A 22 A 23 1 5c = 

A31 A 2 3 ^33 + A33/ 



'fe n + A„ > 

fe 22 + A 22 

; £33 + As3> 



where 5 = 



1 
-A12 



-As 



1 
— A 23 

#33 — £ll &22 — £33 



fell — ^22 £33 — fell 
A 23 



fell — k 2 2 
A 3 1 



fe22 — fe 3 3 

1 



is the transformation 

*' = X (21.18) 



The electro and piezo-optic effects of biaxial, crystals can be handled by 
these infinitesimal transformations, but uniaxial crystals and cubic crystals 
may require finite rotations to re-diagonalize the k~ matrix. In the 

Ik 0\ 

case of cubic crystals we note that I fe A J may be diagonalized by a 

\0 A fe/ 



rotation of 45° about .Vi, giving 



Ik N 

fe + A 
vO fe - A, 



(21.19) 



fe a 2 \ / fe Va? + A'A 



;A 2 Ai fe 



w/aTF a? 
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tan ' — about x 3 and diaonalizes by then rotating through 45° about x 2 
A 2 



70 
giving: 



/k + Va? + A| o > 

k" = k 

\ k - a/A? + Af> 



(21.20) 



The work, can be handled with single column matrices K instead of using 
the square matrices k. If k is a diagonal single column matrix (i.e., the 
single column matrix of ia diagonal matrix), than the almost diagonal 
matrix K + A is diagonalized by the transformation: 



1 








1 


-A* 






1 

A 4 




2A 4 



2Ab 



2Ae 



K\ — Ks K\ — K 2 

— 2A« 

— — 

K 2 - K s Ki- K 2 

-2A4 -2 As Q 

A'o — ^3 ^i — K-z 

-Ae -A 5 



Ki — Kz K 2 — Kz 
A 5 



Ki — Kz 

-A 6 
K x - K 2 K x - K 2 





A 6 



Ac 



K\ — K 2 Ki — Kz 
A* 



K\ — Kz K\ — K 2 

A 5 A4 







K\ — Kz K 2 — Kz 



Ki — Kz 
1 



/ \ 


^1 + A! 


K 2 + A 2 


tfi + Aa 


A 4 


A 5 


A 6 



f^l 


+ Ai 


K 2 


+ A 2 


K a 


+ A 3 


















(21.21) 
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which transforms vectors as x' — 8x where 

A 6 



8 = 



1 
-A B 



A 5 



Ai — A2 Ai — A3 

A4 



K\ — Ki 
-A 5 



1 
-A4 



K.\ — A3 A2 — A3 



A2 — A3 
1 



(21.22) 



If Ai, A 2 and A3 are not all different the preceding analysis falls down as 
some terms become infinite; a finite transformation is needed in this case. 
The difficulty can generally be doged by applying a 45° rotation about one 
of the axes. Sometimes the easiest solution is to rotate through an angle 
<j> about a coordinate axis then solve for the value of <j> that will vanish 
certain terms. As examples of these devices we give the following: 



Kx 
Ki 

Ai 

A., 
A 5 




A, 

A'! 

A 4 







rotated through 
A4 



rotated 45° about X\ becomes 



[K x 


\ 


K x 


+ A 4 


Ai 


- A 4 












1° 


a 



(21.23) 



tan 



£3 is 



about 



Kx 
Ki 
Ki 





Vtf + A£ 




and ro- 
tated 45° 
about xz 

is 



A 


+ Va 4 2 


+ ai 


A 
A 






- Va? 


+ A§ 
























. (21.24) 



D 
E 
k 
K 

a 
a 
X 

€ 

s 
c 

II 



electric induction 
electric field 

dielectric constant matrix (square) 
dielectric constant matrix (single column) 
transformation matrix for vectors 

transformation matrix for tensors (of stress tensor sort) 
stress matrix (single column) 
strain matrix (square) 
strain matrix (single column) 
elastic modulus matrix 
elastic constant matrix 
temperature change of elastic modulus matrix 



11 



BELL SYSTEM TECHNICAL JOURNAL 



h = temperature change of elastic constant matrix 
T s = temperature coefficient of elastic modulus terms 
T c = temperature coefficient of elastic constant terms 
A = temperature expansion matrix 
d — piezo-electric constant matrix 
g = inverse piezo effect matrix 
G = electro-striction matrix 
Z = electro-optic matrix 
it = stress optic matrix 
m — strain optic matrix 

P = polarization = elec. mom. per unit vol. = surface charge per unit 
area 

... P 

7] = susceptibility = -= 



Tabic of Equations 



Transformations 



D = 1-kE 
Air 

e = SX 

X = Ce 

S = C~ l 

h = -C°HC° 

H = -S°hS° 

\l = tAl + PBl 

D = dX 

e = gE 

g = d c 

jT-i = K- 1 ' + ZE 

A'" 1 = A'- 1 " + wX 

A'" 1 = A" 10 + me 

k — J 
P = t—J E= ttE 
4ir 

k - I 



The form * y z 

x' h m\ Hi 

y'/o j»2 n> 

z' h ni3 «3 

h h h 
is the transformation a = m\ m* mz 

Hi «2 «3 

A' = aAa-i 
D' = aD 
E' = aE 
k' = aka c 
K' = a A 
X' = aX 
e' = aj l e 

C = aCcXc 

S' = a^ScT 1 

W = ahcxc 
H = ac l Ha~ l 
ha = CuTCij 
Ha = SljTsa 
d' = ada- 1 
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